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Abstract. In the present contribution, we derive from kinetic theory a unified 
fluid model for multicomponent plasmas by accounting for the electromagnetic 
field influence. We deal with a possible thermal nonequilibrium of the transla- 
tional energy of the particles, neglecting their internal energy and the reactive 
collisions. Given the strong disparity of mass between the electrons and heavy 
particles, such as molecules, atoms, and ions, we conduct a dimensional anal- 
ysis of the Boltzmann equation and introduce a scaling based on the square 
root of the ratio of the electron mass to a characteristic heavy-particle mass. 
We then generalize the Chapman-Enskog method, emphasizing the role of a 
multiscale perturbation parameter on the collisional operator, the streaming 
operator, and the collisional invariants of the Boltzmann equation. The system 
is examined at successive orders of approximation, each of which corresponding 
to a physical time scale. The multicomponent Navier-Stokes regime is reached 
for the heavy particles, which follow a hyperbolic scaling, and is coupled to 
first order drift-diffusion equations for the electrons, which follow a parabolic 
scaling. The transport coefficients are then calculated in terms of bracket op- 
erators whose mathematical structure allows for positivity properties to be 
determined. They exhibit an anisotropic behavior when the magnetic field is 
strong enough. We also give a complete description of the Kolesnikov effect, 
i.e., the crossed contributions to the mass and energy transport fluxes coupling 
the electrons and heavy particles. Finally, the first and second principles of 
thermodynamics are proved to be satisfied by deriving a total energy equation 
and an entropy equation. Moreover, the system of equations is shown to be 
conservative and the purely convective system hyperbolic, thus leading to a 
well defined structure. 



1. Introduction 

Plasmas are ionized gas mixtures, either magnetized or not, that have many 
practical applications. For instance, lightning is a well-known natural plasma and 
has been studied for many years [1] . A second application is encountered in hyper- 
sonic flows; when a spacecraft enters into a planetary atmosphere at hypervelocity, 
the gas temperature and pressure strongly rise through a shock wave, consequently, 
dissociation and ionization of the gas particles occur in the shock layer. Hypersonic 
flow conditions are reproduced in dedicated wind-tunnels such as plasmatrons, arc- 
jet facilities, and shock-tubes [44j [48j E2] ■ A third example was found about two 
decades ago, when large-scale electrical discharges were discovered in the meso- 
sphere and lower ionosphere above large thunderstorms; these plasmas are now 
commonly referred to as sprites 8, 45 1. Fourth, discharges at atmospheric pressure 
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have received renewed attention in recent years due to their ability to enhance the 
reactivity of a variety of gas flows for applications ranging from surface treatment 
to flame stabilization and ignition (see [13J EH EH [53] and references cited therein) . 
Fifth, Hall thrusters are being developed to replace chemical systems for many on- 
orbit propulsion tasks on communications and exploration spacecraft [3J[9j. Finally, 
two important applications of magnetized plasmas are the laboratory thermonuclear 
fusion [6l [49] and the magnetic reconnection phenomenon! in astrophysics |54j . 

Depending on the magnitude of the ratio of the reference particle mean free path 
to the system characteristic length (Knudsen number), two different approaches are 
generally followed to describe the transport of mass, momentum, and energy in a 
plasma [5j : either a particle approach at high values of the Knudsen number (solu- 
tion to the Boltzmann equation using Monte Carlo methods), or a fluid approach 
at low values (solution to macroscopic conservation equations by means of com- 
putational fluid dynamics methods). In this work, we study plasmas that can be 
described by a fluid approach, thus encompassing most of the above-mentioned 
applications. In this case, kinetic theory can be used to obtain the governing con- 
servation equations and transport fluxes of plasmas. Hence, closure of the problem 
is realized at the microscopic level by determining from experimental measurements 
either the potentials of interaction between the gas particles, or the cross-sections. 

A complete model of plasmas shall allow for the following physical phenomena 
to be described 

• Thermal non equilibrium of the translational energy, 

• Influence of the electromagnetic field, 

• Occurrence of reactive collisions, 

• Excitation of internal degrees of freedom. 

So far, no such unified model has been derived by means of kinetic theory. Besides, a 
derivation of the mathematical structure of the conservation equations also appears 
to be crucial in the design of the associated numerical methods. Based on our 
previous work, we investigate in the present study the thermal noncquilibrium of the 
translational energy [39] and the influence of the magnetic field [27] ■ We generalize 
the Chapman-Enskog method within the context of a dimensional analysis of the 
Boltzmann equation, emphasizing the role of a multiscale perturbation parameter 
on the collisional operator, the streaming operator, and the collisional invariants of 
the Boltzmann equation. Then, we obtain macroscopic equations eventually leading 
to a sound entropy structure. Moreover, the system of equations is shown to be 
conservative and the purely convective system hyperbolic. Let us now describe in 
more detail how these issues are currently addressed in the literature. 

First, a multiscale analysis is essential to resolve the Boltzmann equation govern- 
ing the velocity distribution functions. We recall that a fluid can be described in the 
continuum limit provided that the Knudsen number is small. Besides, in the case 
of plasmas, a thermal nonequilibrium may occur between the velocity distribution 
functions of the electrons and heavy particles (atoms, molecules, and ions), given 
the strong disparity of mass between both types of species. Ergo, the square root of 
the ratio of the electron mass to a characteristic heavy-particle mass represents an 
additional small parameter to be accounted for in the derivation of an asymptotic 
solution to the Boltzmann equation. Literature abounds with expressions of the 
scaling for the perturbative solution method. For instance, significant results are 
given in references [TBI [TBI E21 EH [55] . Yet, Petit and Darrozes [46] have suggested 
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that the only sound scaling is obtained by means of a dimensional analysis of the 
Boltzmann equation. Subsequently, Degond and Lucquin pjJl [20] have established 
a formal theory of epochal relaxation based on such a scaling. In their study, the 
mean velocity of the electrons is shown to vanish in an inertial frame. Moreover, the 
heavy-particle diffusive fluxes were scarcely dealt with since their work is restricted 
to a single type of heavy particles, and thus no multicomponent diffusion was to be 
found; in such a simplified context, the details of the interaction between the heavy 
particles and electrons degenerate and the positivity of the entropy production is 
straightforward. We will establish a theory based on a multiscale analysis for a 
multicomponent plasma (which includes the single heavy-particle case) where the 
mean electron velocity is the mean heavy-particle velocity in absence of external 
forces. As an alternative, Magin and Degrez [39] have also proposed a model for a 
multicomponent plasma based on a hydrodynamic referential. They have applied 
a multiscale analysis to the derivation of the multicomponent transport fluxes and 
coefficients. However, the proposed treatment of the collision operators is heuris- 
tic. Moreover, since the hydrodynamic velocity is used to define the referential 
instead of the mean heavy-particle velocity, the Chapman-Enskog method requires 
a transfer of lower order terms in the integral equation for the electron perturbation 
function to ensure mass conservation. Finally, we also emphasize that the develop- 
ment of models for plasmas in thermal equilibrium shall always be obtained as a 
particular case of the general theory. 

Second, the magnetic field induces anisotropic transport fluxes when the elec- 
tron collision frequency is lower than the electron cyclotron frequency of gyration 
around the magnetic lines. Braginskii [lOj has studied the case of fully ionized 
plasmas composed of one single ion species. Recently, Bobrova et al. have gener- 
alized the previous work to multicomponent plasmas. However, the scaling used 
in both contributions does not comply with a dimensional analysis of the Boltz- 
mann equation. Lucquin [36|, 137] has investigated magnetized plasmas in the latter 
framework. Nevertheless, the same limitation is found for the diffusive fluxes as in 
reference [THl HO] • Finally, Giovangigli and Graille [27] have studied the Enskog ex- 
pansion of magnetized plasmas and obtained macroscopic equations together with 
expressions of transport fluxes and coefficients. Unfortunately, the difference of 
mass between the electrons and heavy particles is not accounted for in their work. 

Third, plasmas are strongly reactive gas mixtures. The kinetic mechanism com- 
prises numerous reactions [12] : dissociation of molecules by electron and heavy- 
particle impact, three body recombination, ionization by electron and heavy-particle 
impact, associative ionization, dissociative recombination, radical reactions, charge 
exchange. . . Giovangigli and Massot [28] have derived a formal theory of chemi- 
cally reacting flows for the case of neutral gases. Subsequently, Giovangigli and 
Graille [27] have studied the case of ionized gases. We recall that their scaling 
does not take into account the mass disparity between electrons and heavy parti- 
cles. Besides, in chemical equilibrium situations, a long-standing theoretical debate 
in the literature deals with nonuniqucness of the two-temperature Saha equation. 
Recently, Giordano and Capitclli [29 have emphasized the importance of the phys- 
ical constraints imposed on the system by using a thermodynamic approach. A 
derivation based on kinetic theory should further improve the understanding of the 
problem. Choquet and Lucquin [15] have already studied the case of ionization 
reactions by electron impact. 



4 BENJAMIN GRAILLE, THIERRY E. MAGIN*, AND MARC MASSOT 

Fourth, molecules rotate and vibrate, and moreover, the electronic energy lev- 
els of atoms and molecules are excited. Generally, the rotational energy mode 
is considered to be fully excited above a few Kelvins. In a plasma environment, 
the vibrational and electronic energy modes are also significantly excited. The de- 
tailed treatment of the internal degrees of freedom is however beyond the scope of 
the present contribution where we will only tackle the translational energy in the 
context of thermal nonequilibrium. The reader is thus referred to the specialized 
literature [TTJ IM EU • 

Fifth, the development of numerical methods to solve conservation equations re- 
lies on the identification of their intrinsic mathematical structure. For instance, the 
system of conservation equations of mass, momentum, and energy is found to be 
nonconservative for thermal nonequilibrium ionized gases. Therefore, this formula- 
tion is not suitable for numerical approximations of discontinuous solutions. Coqucl 
and Marmignon [17] have derived a well-posed conservative formulation based on a 
phenomenological approach. Nevertheless, their derivation is not consistent with a 
scaling issued from a dimensional analysis. We will show that kinetic theory, based 
on first principles, naturally allows for an adequate mathematical structure to be 
obtained, as opposed to the phenomenological approach. 

In this work, we propose to derive the multicomponcnt plasma conservation 
equations of mass, momentum, and energy, as well as the expressions of the as- 
sociated multicomponent transport fluxes and coefficients. The multicomponcnt 
Navier-Stokes regime is reached for the heavy particles, which follow a hyperbolic 
scaling, and is coupled to first order drift-diffusion equations for the electrons, which 
follow a parabolic scaling. We deal here with first-order equations, thus one order 
beyond the expansion commonly investigated for the electrons in the literature. 
The derivation relies on kinetic theory and is based on the ansatz that the particles 
constitutive of the plasma are inert and only possess translational degrees of free- 
dom. The electromagnetic field influence is accounted for. In Section [21 we express 
the Boltzmann equation in a noninertial reference frame. We show that the mean 
heavy-particle velocity is a suitable choice for the referential velocity. This step is 
essential to establish a formalism where the electrons follow the bulk movement of 
the plasma. Then, we define the reference quantities of the system in order to derive 
the scaling of the Boltzmann equation from a dimensional analysis. The multiscale 
aspect occurs in both the streaming operator and collision operator of the Boltz- 
mann equation. Consequently, Section [3] is devoted to the scaling of the partial 
collision operators between unlike particles. Besides, we determine the space of 
collisional invariants associated with respectively the electrons and the heavy par- 
ticles. In Section IU we resort to a Chapman-Enskog method to derive macroscopic 
conservation equations. The system is examined at successive orders of approxima- 
tion, each of which corresponding to a physical time scale. For that purpose, scalar 
products and linearized collision operators are introduced. The global expressions 
of the macroscopic fluid equations are then provided up to Navier-Stokes equations 
for the heavy particles and first-order drift-diffusion equations for the electrons. We 
also prove that our choice of referential is essential in order to reach this expansion 
level. In Section [H we establish the formal existence and uniqueness of a solution 
to the Boltzmann equation. The multicomponent transport coefficients are then 
calculated in terms of bracket operators whose mathematical structure allows for 
the sign of the transport coefficients to be determined; in particular, the Kolesnikov 
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effect, or the crossed contributions to the mass and energy transport fluxes coupling 
the electrons and heavy particles. The explicit expressions of the transport coef- 
ficients can be obtained by means of a Galerkin spectral method [14] disregarded 
in the present contribution. Finally in Section [6] the first and second principles 
of thermodynamics are proved to be satisfied by deriving a total energy equation 
and an entropy equation. Then, we establish, from a fluid standpoint, a conser- 
vative formulation of the system of macroscopic equations. We also identify the 
mathematical structure of the purely convective system. Hence, we demonstrate 
that kinetic theory shall be used as a guideline in the derivation of the macroscopic 
conservation equations as well as in the design of the associated numerical methods. 

Beyond the obvious interest of such a study from the point of view of the ap- 
plications and design of numerical schemes, the present contribution also yields a 
formal kinetic theory of mixtures of separate masses, where the light species obey a 
parabolic scaling whereas the heavy species obey a hyperbolic scaling. The original 
treatment of the two different scalings for fluid flows was first provided by Bardos 
el al. 2J. In their study, the purely hyperbolic scaling yields the compressible gas 
dynamics equations, whereas the purely parabolic scaling leads to the low Mach 
number limit. These scalings are quite classical and both of them can be used for 
various asymptotics such as the Vlasov-Navier-Stokes equations in different regimes 
investigated by Goudon et al. [32(32]. Yet, a rigourous derivation of a set of macro- 
scopic equations in the situation where the hyperbolic and parabolic scalings are 
entangled in the same problem is an original result obtained in the present work. 

2. BOLTZMANN EQUATION 

2.1. Assumptions. 

(1) Our plasma is described by the kinetic theory of gases based on classical 
mechanics, provided that: a) The mean distance between the gas particles 
l/(n ) 1 / 3 is larger than the thermal de Broglie wavelength, where n° is a 
reference number density, b) The square of the ratio of the electron thermal 
speed Vg to the speed of light is small. 

(2) The reactive collisions and particle internal energy arc not accounted for. 

(3) The particle interactions are modeled as binary encounters by means of 
a Boltzmann collision operator, provided that: a) The gas is sufficiently 
dilute, i.e., the mean distance between the gas particles l/(n ) 1 / 3 is larger 
than the particle interaction distance (cr ) 1 / 2 , where a is a reference dif- 
ferential cross-section common to all species, b) The plasma parameter, 
quantity proportional to the number of electrons in a sphere of radius equal 
to the Debye length, is supposed to be large. Hence, multiple charged par- 
ticle interactions are treated as equivalent binary collisions by means of a 
Coulomb potential screened at the Debye length [1] [22] . 

(4) A plasma is composed of electrons and a multicomponent mixture of heavy 
particles (atoms, molecules, and ions). The ratio of the electron mass m® 
to a characteristic heavy-particle mass m° is such that the nondimensional 
number e = */ mj? /m? is small. 

(5) The number of Mach, defined as a reference hydrodynamic velocity divided 
by the heavy-particle thermal speed Mh = v° /V®, is supposed to be of the 
order of one. 
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(6) The macroscopic time scale t° is assumed to be comparable with the heavy- 
particle kinetic time scale i° divided by e. The macroscopic length scale is 
based on a reference convective length L° — v°t°. 

(7) The reference electrical and thermal energies of the system are of the same 
order of magnitude. 

The mean free path 1° and macroscopic length scale L° allow for the Knudsen 
number to be defined Kn = l°/L°. It will be shown that this quantity is small, 
provided that assumptions (JU)-© are satisfied. Therefore, a continuous description 
of the system is deemed to be possible. 

2.2. Inertial frame. The choice of an adequate referential will prove to be essential 
in the following multiscale analysis. Two referential are commonly used in the 
literature. Degond and Lucquin [101 120] work in the inertial frame, as Ferziger and 
Kaper [25] . The second referential is presented in the following section. Considering 
assumptions CO)-©, the temporal evolution of the velocity distribution function f* 
of the plasma particles i is governed in the phase space (x*, c*) by the Boltzmann 
equation [T3"I[2"5] 

(2.1) W)=#, *eS, 

where symbol S is the set of indices of the gas species. Dimensional quantities are 
denoted by the superscript *. The streaming operator reads 

(2.2) @*(f*) = dt .f* + c l d ^f* + 4L(E* + c*AB*yd ct f*, *6S, 

in an inertial frame. Symbol t* stands for time, E*, the electric field, B*, the 
magnetic field, m*, the mass of the particle i, and q*, its charge. The collision 
operator is given by 

(2.3) # = XX (//,//), ieS, 

with the partial collision operator of particle j impacting on particle i 

(2.4) a*, (f*, //) = J (f*'f*> - /,*//) K - 4l4 dwdc i> i,j e s. 

After collision, quantities are denoted by the superscript '. The differential cross- 
section <7y —<^ij [n*j\ c i — C ^\ 2 /(^-bT ), w-e] depends on the relative kinetic energy 
of the colliding particles and the cosine of the angle between the unit vectors of 
relative velocities u; = (c*' — c*')/\c*' — c*'| and e = (c* — c*)/\c* — c*\. Quantity 
fi*j = m*m*/ (m*+m*) is the reduced mass of the particle pair, T , a reference tem- 
perature, and ke, Boltzmann's constant. Therefore, the differential cross-sections 
are symmetric with respect to their indices i,j € S, i.e., a*j = a* { . 

2.3. Noninertial frame. Sutton and Sherman |51j. as Chapman and Cowling [14j. 
have proposed a noninertial frame based on the hydrodynamic velocity 

(2.5) p*v* = W m*c*f*dc*, 

jes J 

where the mixture mass density is defined as p* — YljesPj' Symbol p* — n*m* 
stands for the partial mass density, and n* — J f*dc*, the partial number density. 
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It is a convenient choice since it is the referential associated with the definition of 
the peculiar velocities 

(2.6) Cr =c*i-v*, ie S, 

induced from the usual momentum constraint. We infer from definition (|2.5[) that 
the global diffusion flux vanishes 



(2.7) ^ / m;e;'« ./;.;,;; = o, 

jes J 

that is, the standard momentum constraint. 

Given the strong disparity of mass between the electrons and heavy particles, 
a frame linked with the heavy particles appears to be a rather natural choice for 
plasmas, as fully justified in the following detailed multiscale analysis. Thus, we 
define the mean electron velocity and mean heavy-particle velocity 

(2.8) px = / K<f e *dc*, P %v* h = W »';<./; '!<•;• 

where the heavy-particle mass density reads p* h — J^jeaPj- ^ n ^ nc v h referential, 
the free electrons interact with heavy particles whose global movement is frozen 
in space. A similar viewpoint is commonly adopted in the quantum theory of 
molecules when the Born-Oppcnhcimcr approximation is used to study the motion 
of the bound electrons about the nuclei [7]. Based on the following definition of 
peculiar velocities 

(2.9) C* = c$-v* h , zeS, 
the heavy-particle diffusion flux is shown to vanish 



(2.io) I <»: c ;i:<k 



0. 

( J J "J J 

For now, we defer the choice of the referential velocity. Therefore, we use the 
symbol u* to define the peculiar velocities C™* = c* — it*, i G S. Then, the 
Boltzmann equation is expressed in a frame moving at u* velocity by means of the 
latter change of variables. Hence, the streaming operator (|2.2p is transformed into 
the expression 

(2.11) 

mur) = d*ft + (cr + u*) -d x *ft + 4 [e* + (cr + «*) ^ -d cr ft 

m i 
Dn* 

where D/Di* = df + u* d x *. The partial collision operator (|2.4[) is found to be 

(2.12) Sti U*= #) = / (r //* - fiff) \Cr exclude?*, i,j e S. 

In a noninertial frame, the velocity distribution function /*, the differential cross- 
section a*j = a*j [/j,*j\C™* - C^| 2 /(kBr ), o>-e], as well as both the unit vectors 
w = (C? 1 * -C?'*)/\C?'* -C?'*\ and e = {C?* -C?*)/\C?* -C?*\ depend on the 
peculiar velocities. Nevertheless, for simplicity reasons, notations are unchanged 
with respect to the inertial frame, where the previous quantities depend on the 
absolute velocities. Moreover, we define collisional invariants. 
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Definition 2.1. The space of scalar collisional invariants I"* is spanned by the 
following elements 



(2.13) 



< V u '" S+l/i = ("W) ieS , ^£{1,2,3}, 

where symbol rfi denotes the cardinality of the set of species S. 
Besides, a scalar product is introduced 

(2.i4) «e\or = E / ^©c/ dc ; 



1 11 A 



for families £* = (Cf)igs and = (C^)ieS - Symbol stands for the maximum 
contracted product in space and symbol ~ for the conjugate transpose operation. 
Ergo, the collision operator 3* — (3*) ie s defined in eq. (|2.3p obeys the following 
property. 

Property 2.1. The collision operator 3* is orthogonal to the space of collisional 
invariants I u * , i.e., ({i(> u > 1 *, 3*})™* = 0, for all I G {1, . . . ,n s +4}. 

Proof. The projection of the collision operator 3* onto ifj u ' 1 * , i e {1,... , tt s +4}, is 
shown to be 

see for instance Chapman and Cowling [14] . The latter expression vanishes for all 
I £ {l,...,n s +4}. □ 

Finally, the macroscopic properties can be expressed by means of the scalar 
product of the distribution functions and the collisional invariants 

f <(/*,#*»* = pi ieS, 

«r>« s +-r = p*(v* v -u*), ^e{i,2,3}, 

i ((/*,^" S+4 *F = p*e* + \p*{v* - v.*)- (v*-u*), 
where quantity e* stands for the gas thermal energy per unit mass. 

2.4. Dimensional analysis. A sound scaling of the Boltzmann equation is de- 
duced from a dimensional analysis inspired from Petit and Darrozes |46j . First, 
reference quantities are introduced in Table [T] The characteristic temperature, 
number density, differential cross-section, mean free path, macroscopic time scale, 
hydrodynamic velocity, macroscopic length, and electric and magnetic fields are 
assumed to be common to all species. The nondimensional number 

(2.15) e = 

quantifies the ratio of the electron mass to a reference heavy-particle mass. Ac- 
cording to assumption ((U, the value of e is small. Consequently, electrons exhibit 
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Common to all species 

Temperature T° 

Number density n° 

Differential cross-section a 

Mean free path 1° 

Macroscopic time scale t° 

Hydrodynamic velocity v° 
Macroscopic length 

Electric field E° 

Magnetic field B° 

Electrons Heavy particles 

Mass ml? m° h 

Thermal speed V° 14° 

Kinetic time scale t° 



a larger thermal speed than that of heavy particles 



Moreover, the electron and heavy particle temperatures may be distinct, provided 
that eq. (|2.16|) does not fail to describe the reference thermal speeds. The differential 
cross-sections are of the same order of magnitude a . Hence, the characteristic mean 
free path 1° = l/(n°a°) is found to be identical for all species. As a result, the 
kinetic time scale, or relaxation time of a distribution function towards its respective 
quasi-equilibrium state, is lower for electrons than for heavy particles 

rO 70 f o 

4 = ^ = |- 

Assumption $Q states that the macroscopic time scale reads 

t° 

(2.18) t° = A 

e 

It is shown in Section 2] that this quantity corresponds to the averaged time during 
which electrons and heavy particles exchange their energy through encounters. In 
addition, the macroscopic temporal and spatial scales are linked by the expression 

(2.19) L° = v°t a , 

where the hydrodynamic velocity is determined by the Mach number My t — v /Vff. 
Given assumption (O, the Mach number is of the order of one. Hence, the Knudsen 
number 

(2.20) Kn = T° = -k- 

is small. Finally, following assumption (J7|, the reference electric field is such that 

(2.21) q°E°L = k B T°. 
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The intensity of the magnetic field is governed by the Hall numbers of the electrons 
and heavy particles 

o R o o R o 

(2.22) A= 9_f_ t o =e i-6 j ^=«4-t| \=efi , 

m% rrvi 

defined as the Larmor frequencies, q°B° /m e for the electrons and q°B° /mj for the 
heavy particles, multiplied by their corresponding kinetic time scale. The magnetic 
field is assumed to be proportional to a power of e by means of an integer b < 1. 
The physical interpretation of the b parameter appears at the end of Section [5] 

The dimensional analysis can be summarized as follows: a) Two spatial scales 
were introduced, one spatial scale at the microscopic level and one spatial scale 
at the macroscopic level; b) Whereas three temporal scales were defined, two time 
scales at the microscopic level, respectively for the electrons and for the heavy 
particles, and one time scale at the macroscopic level, common to all species. 

Nondimcnsional variables are based on the reference quantities. They are de- 
noted by dropping the superscript *. In particular, one has the following expressions 
for the particle velocities 

(2.23) < = V c °c c , c$ = V£a, i £ H, 

where symbol H stands for the set of indices of heavy particles. Both the reference 
hydrodynamic velocity and mean heavy-particle velocity are equal to v°. Indeed, 
the hydrodynamic velocity defined in eq. (|2.5p is found to be 

(2.24) (p h + e 2 p c )M h v = p h M h v h + s 2 p c v c , 

whereas the mean electron and heavy-particle velocities given in eq. (|2.8|) read 

(2.25) p c v c = i J c c f c dc c , p h M h v h J rrijCjfjdcj. 

The peculiar velocities are given by the relations 

(2.26) C™ = c c eM h u, C™ = C* - M h u, i e H. 

Usually, they are associated with the momentum constraints of the mixture, so that 
the natural choice is u = v. In such a case, we get the following relation 



E 



J mjCJfjACJ+e J C* '£dC£ = 0. 



However, the hydrodynamic velocity of the mixture, electrons included, can also 
be expanded in the e parameter and thus receives contributions at various e orders 
in the Chapman-Enskog method. Since the change of referential should not differ 
depending on the expansion order, we could mimic the approach of Lucquin and 
Degond [19, 20, 36 and take u — 0, which means working in the incrtial framework. 
However, we follow a different path not only by choosing the mean heavy-particle 
velocity as referential velocity, u — Vh, but also by defining the peculiar velocities 
based on this quantity, as opposed to Petit and Darrozes 46J. The rationale for 
such a choice is threefold: a) This quantity does not depend on e while still being 
a perturbation of the hydrodynamic velocity of the complete mixture up to second 
order in e 

(2.27) (p h + e 2 Pc )M h (v - v h ) = e J f c AC^ , 
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since quantity / C™ h f e &C™ h taken on a perturbation of a Maxwell-Boltzmann dis- 
tribution will be of 0(e) in the framework of the Chapman-Enskog expansion pre- 
sented in Section 31 b) It will prove to be the natural referential in which the heavy 
particles thermalize in the context of the proposed multiscale analysis; c) It will 
also prove to be the only available choice for electron thermalization and successive 
order resolubility, thus making the proposed change of referential optimal and lead- 
ing to a rigourous framework as well as a simplified algebra. In the sequel, since 
there is no ambiguity, we will drop the Vh superscript in the use of the peculiar 
velocities C" h and C/ h , i G H. 

Consequently, the heavy-particle diffusion flux vanishes, as shown in eq. (|2.10|) 

(2.28) I m iQfj dC j =°- 

Thus, the Boltzmann equation (|2.1[) can be expressed in nondimensional form, 
respectively for the electrons and heavy particles, as 



(2.29) d t f c + j^(C c + eM h v h )-d x f c + e- (1+b \ e [(C c + eM h v h )^B] -d c J c 



G 



1 q e E - sM h ^) -del - d c J c ®C c :d x v h = ±3 C 



(2.30) d t fi + ^(d + M h v h )-d x f t + e 1 -"^-^ + M h v h )^B] -d c Ji 

+ (t^^tE - M h ^f) -dej - dcJivQ-.dvVh = i G H, 



where the collision operators read 

(2.31) 3 C = 3 oc 3c, (/c Jj) , 

(2.32) 3i = |afe(JS,JL) + X;^U./i). ieH - 

The collisional invariants (|2. 13|) depend on the mass ratio as well, as shown in their 
nondimensional form. 

Definition 2.2. The space of scalar collisional invariants X is spanned by the 
following elements tp l £ = (tp l e , tp l h ), I G {1, . . . , 77 s +4}, with 



(2.33) 



It is worth noticing the influence of the hierarchy of scales; whereas the scaling 
does not introduce any structural change in the mass and energy collisional invari- 
ants, the electron contribution disappears from the momentum collisional invariant 
vector in the limit of e tends to zero. A similar behavior can be observed for the 
total mass; however, the single species collisional invariants are not affected. 

Let us underline that eq. (|2.29[) for the light species is typical of a parabolic scal- 
ing, which corresponds to the low Mach number limit for the electron gas, whereas 
eq. (|2.30[) for the heavy species is typical of a hyperbolic scaling, which corresponds 
to the compressible gas dynamics for the heavy-species gas mixture [2] . The present 





e 2 5 ■ 


< = 


(mi%) ieH , 




c s+ " = 








v G {1,2,3} 


c s+4 = 


2 1 


= 


(>;a-a) ieH . 
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scaling is thus intermediate between the usual cases and the mathematical structure 
of the resulting system of macroscopic equations has to be identified. 

For a family £ = (£i) igS , we introduce two separate contributions: £ e , concerning 
the electrons, and £/j = (£i) i6H , concerning the heavy particles. Hence, the scalar 
product between the families £ = (£j)ies an d C = (Ci)ieS defined in eq. (|2.14p is 
decomposed into a sum of partial scalar products with different scales 

(2.34) «e,c» = «ee,a +e 3 <(&»Cfc)>h. 

given by the expressions 

(2.35) Ce)> a = / CeOCo dC c , = E / 0©6 dQ. 

Finally, we introduce the collision operator 3 e = (e<Jei where eq. (|2.29[) has 

been multiplied by a factor e 3 corresponding to a coherent scaling of the two Boltz- 
mann equations. Then, we derive the following property. 

Property 2.2. The collision operator 3 S is orthogonal to the space of collisional 
invariants I, i.e., ((ip l £ ,3e)) = 0, for all I 6 {f , . . . ,?t s +4}. Furthermore, the terms 
of ((ip e ,3e)) cancel by pair of interaction, i.e., 



(2-36) «^Uee»e = 0> 

(2-37) £<^Uei>> e + «C0 fc e)) h =O, 

(2-38) J ' eH ^2U l h ,3 hj )) h = 0, 



respectively for the electron, electron heavy-particle, and heavy-particle interactions. 

Proof. The projection of the collision operator 3e onto tfj l e , I £ {1, . . . , n s +4}, is 
given by the expression 

((4,3 e )) = eU l c , 3 CC )) C + e £ Ul 0ei» e + eUl 3h*)) h + e 2 £ «^U*i» h - 

jGH jGH 

The terms of this sum are examined by interaction pairs 

^)|C -C cl | CTccl du;dC c dC ol , 

^OICe-eQIa^dwdCedC., 



^■OlQ-Ql^dwdadCj. 

These expressions vanish and thus the sum ((ip l e ,3s)} =0. □ 

The multiscale analysis occurs at three levels: a) In the kinetic equations (|2.29|) 
and (|2.30| : b) In the collisional invariants (|2.33|) and thus in the conservation 
of the associated macroscopic quantities; c) In the collision operators. Actually, 
encounters between particles of the same type are dealt with as usual, whereas the 



«Vi 3cc» c = \ J (/ e 7el - fc UM + 4l ~ $ ~ 

= ^£ hfifj-ufi)M+$-<- 

£ #U«»„ 
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collision operators between unlike particles depend on the e parameter via their 
relative kinetic energy and velocity, and the vectors u) and e. The scaling of these 
operators is investigated in the following section. 

3. Preliminary results 

3.1. Electron heavy-particle collision dynamics. The study of the electron 
heavy-particle collision dynamics yields the dependence of the peculiar velocities 
on the e parameter. First, we express the momentum conservation in terms of 
the peculiar velocities. Considering a collision of a heavy species, i £ H, against 
an electron, the peculiar velocities after collision C[ and C e are related to their 
counterpart before collision C$ and C c 



(3.1) 



F 777 " f 

C[= ~^C C + ^Ci+s —\eCi-Q>\u, ieH, 

rrii + e z m,- + e rrii + £ 

C n = ■ — - C r . H ■ — - d -s e C,; -C t w, 



mj + £ z mj + £ rrii + e z 

where the direction of the relative velocities after collision is defined in their center 
of mass by 

eC! - CL 



\eQ - C c l 

Symbol s stands for an integer either equal to +1 for the collision operator die, 
i £ H, or — 1 for 3d, i S H. This notation allows for the interaction considered 
in eq. (|2.12|) to be distinguished. We are now able to expand the crossed-collision 
operators. 

3.2. Expansion of the collision operator die- The dimensional analysis yields 
the following expression of the nondimensional collision operator Sic, i £ H, 

(3.2) 3ie(fi,fe)(C i ) = J a le (| 7c | 2 ,u;.^) | e Q - C e \ 

"j5 (C/)/ c (CO - fi (Q)/ e (C e )] do, dC c , 

where the relative kinetic energy and the vector e are expressed by means of the 
vector 7 = s(eC l - C e )/(1 + e 2 /m i ) 1 / 2 . 

We then introduce the generalized momentum cross-section [14] in a thermal 
nonequilibrium context 

PIT 

(3.3) QtHhd 2 ) = 2tt / a lo (| 7o | 2 ,cos0)(l - cos' 0) sin0 d0, ieH, 

Jo 

where symbol 9 stands for the angle between the vectors u> and e. For I = 1, this 
cross-section represents the average momentum transfered in encounters between i 
heavy particles and electrons for a given value of the relative kinetic energy. 

Theorem 3.1. The collision operator Sic, ieH, can be expanded in the form 

(3.4) 3 l c(f t J c KC l ) = e8i(f l Jc)(C t )+e 2 Sl(f l Jc)(C t )+e 3 SUf l Jc)(C t )+0(e 4 ). 

The zero-order collision operator S%{fi , fc){Ci), ieH, vanishes. The first- order 
term S} e , i e H, reads 

(3-5) &U >£)(«) —dcji (Ci)- f Q\ e \ft e \ 2 ) |7cl7e /e(7e) d7c * 6 H. 
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The second order term 3 2 e , i 6 H, is found to be 

(3.6) aiU,/ e )(Q) = (<*)<*): / Q^{\lf)\l,\dcMl e )^ e d 7c 

+ ^ jG / i (C i ): | 0L 2) (|7 c | 2 )l7cl(l7 c | 2 I-3 7 e«7c)/c(7o)d7o 

+ -^ Q .£(G): / gL 1) (|7c| 2 )l7cl7e^7e/e(7e)d7c- 
Finally, the third-order term 3f e; z G H, is given by 

(3-7) = ^dcAUCiVCi)® f QilHhf) |7el^ e C e /e(7e)®7e d 7o 

+ ^ a / i (a> | gL 1) (l7c| 2 )l7el7c«7c-^/c(7e)d7 c 
+ ^dcAdcJiiCi)®^)® I 0L 1) (|7c| 2 )l7el7c«7c^c o /c(7c)d7 c 
+ ^dc t {dcMC % )®C t )Q J QL 2) (|7c| 2 )l7cl(l7c| 2 I-37c«7c)^ Q ./c(7c) d 7o 

-i^aaiW)© J gL 1) (l7e| 2 )l7el(l7e| 2 I + 7e«7e)«7e/e(7c)d7e 

-^ ac ,/,(a)0 1 gL 2) (|7c| 2 )l7cl(l7c| 2 I-37c^7c)«7c/c(7e)d7c 

+ I ^2^ iCiCi / i (a)0 | gL 3) (|7e| 2 )l7el(3|7c| 2 I-57e^7c)^7e/c(7c)d7e 

+ ^-3L(/ I Jc)(C 1 ). 

Proof. The change of variable dC = — (1 + e 2 /mi) 3 / 2 d 7c allows for the differential 
cross-section dependance on e to be eliminated 



M/*,/e)(G) = | ^J 2 ,^) |7el(l + £ 2 / 



2 

TO,;) 



/ i (c;)/ c (c c ')-./ 4 (a)/e(c e ) 



du> d7 e , ieH. 



Then, the peculiar velocities are expanded in a power series of e 

C/ = a+£^-a 1 -e 3 ^ ? a 1 + 0(e 4 ), a x = - 7c + | 7 >, ieH, 
C£ = -| 7e |w + ea + e 2 ^-a 2 + 0(s 4 ), a 2 = - 7e + || 7 >, 
Ce = - 7c + ea-e 2 ^- 7o + 0(e 4 ). 

Hence, the distribution functions are found to be 

/, {CD = f (d) + e^-dcji (CO-oi + e 2 ^d 2 CiCi fi (Ci):oi®oi 

+ ^e^^QC,/. (C,)0ai®ai®ai - e 3 ^d c J t (Q)-ai + C(e 4 ), i e H, 
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fo (CO = / c (-| 7e |«) + edcj c (-|7 e |w)-C4 + e 2 \d 2 CcC J c (-|7.M:«»a 

+ ^^-^ cCe / (-|7 c |^):C^a 2 + 0( £ 4 ), 
/c (C„) = / e (-7c) + ed c J c (-7 C )-Q + £ 2 i^ cCo / c (-7 e ):Ci®Q 

- e2 2^- 9 c /c (-7c)7c + £% CeCe / e (-7 e )©a®a®a 

- ^ ^hcje (-7c):a^7e + 0(£ 4 )- 

Combining these equations, the zero-order term J? e , i G H, is thus given by 
a?e(i J /e)(a) = / i (Q) | a le (| 7c | 2 ,^e) 

l7o| 3 [/e ("|7el«) - /e HTek)] de d| 7o |. 

Intertwining e with u;, the integral is shown to vanish. Then, eqs. (j3.5|) - (|3.7() are 
obtained after some lengthy calculation. □ 

Theorem 13.11 admits three corollaries. 

Corollary 3.1. The first-order collision operator d} e (fi,fe)i i£H, vanishes when 
f e is an isotropic function of the velocity C c . 

Proof. Expression (|3.5|) immediately yields that the integrand is an odd function of 
7 e if f e is isotropic in the mean heavy-particle frame, so that the first-order collision 
operator vanishes. □ 

Remark 3.1. So far, we note that this property is strongly related to our choice 
of referential. For example, such a property is not satisfied when u = 0. Thus, the 
structure of the expansion of collisional integrals depends on the initial choice of 
referential. We will come back to this point in Section 14751 

A collision frequency is defined as a Maxwell-Boltzmann averaged momentum 
cross-section 



^/QL 1) (l7e| 2 )|7e| 3 /e°(7e)d7 e! i G H, 



T c _ 

where / c °( 7o ) = n c exp h 7o 7 /(2T )] / (2irT c f 2 . 

Corollary 3.2. If f c ° = n c cxp [-C C -C C /(2T C )] / (2ttT c ) 3/2 , the second-order colli- 
sion operator reads 

(3-8) 8 2 c (f l J c °KC l ) = ^(d Ci if l C t ) + ^A c jX i e H. 

6mi \ rrii / 

Proof. A direct calculation of 3 2 c (fi , f c °)(Ci) given in (|3.6p immediately yields ex- 
pression (JMJ) if / c ° = n c cxp [-C e -C e /(2T e )] / (2ttT c ) 3/2 . □ 

Corollary 3.3. The third-order collision operator S 3 c (fi , f e )> i G H, vanishes when 
f e is an isotropic function of the velocity C c . 

Proof. Expression (|3.7[) immediately yields that the integrand is an odd function 
of 7 C if / c is isotropic in the mean heavy-particle frame, so that the third-order 
collision operator vanishes. □ 
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3.3. Expansion of the collision operator 3d- The dimensional analysis yields 
the following expression of the nondimcnsional collision operator 3d, i S H, 



(3.9) dei(fe,fi)(C e 



. . — )|C e eCi\ 

feiCiMiCD-feiaUiCi) 



do; dd. 



The original set of variables {C c , Ci,uj} is retained. We introduce the momentum 
cross-section 

(3.10) Qe^dCel 2 ) = 2tt / a cl (\C c \ 2 ,cos6)(l-cos9)sm9d9, igH, 

Jo 

representing the average momentum transfered in encounters between electrons and 
heavy particles ieH. It is equal to the cross-section Q^} ■ 

Theorem 3.2. The collision operator 3d, i£H, can be expanded in the form 

(3.11) 3ei(fe , /0(C e ) = 3° ct (f c , /0(C e ) + £^(/c , /0(C e ) + , fi)(C e ) 

+ s"3Uf c ,.mC c ) + 0(e 4 ). 



The zero-order term 3%' * G H, is given by the expression 
(3.12) 



3° ei (fe,fi)(C e )= J A(G)dCi J a CJ (|C c | 2 , 

The first-order term 3] c , ieH, reads 



f c (\C c \u>)-f e (C e ) 



(3.i3) 4(/ e , fi )(c e ) = (| /i (a)a da) • 

|a Co | a ej (|C c | 2 ,|^-a;)[/ c (C c )-/ c (|C c |u;)]|C c |du; 

+ | ^(ICel 2 , gy-wJICel [d Ce / e (|C7 |a») - 9 Ce / e (C e )] dwj. 

77ie second order term 3 2 e , i £ H, is found to be 
(3.14) 

3e 2 »(/e , /»)(Ce) = ^^(G.) | /, (Q) dQ + (C e ) : | J? (a)C^Q dQ, 



with 



K 2 C :\C C ) = d Ce - J a ct (\C c \ 2 ,^)(C c |C e |a,)|C e |/ e (|C e |a,) du 

-\C e \C e - J d Ca a cl (\C c \ 2 ,^-u;)[f c (\C c \u>) - f c (C c )] du>, 
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K 2 c f(C c ) = dh eCe J a cl (\C c \ 2 , [/ e (Ct|«) - / e (C e )] da; 

+ 2 /" a Cc (a cl (|C c | 2 ,|^-a;)|C |)®[a c J c (C )-ac c /c(|C |a;)] da; 

+ |C | | Miar, 1 g r -w)[^ ct A(a ) )-^. ai A(ia,|w)] da; 

+ 2|C C | /" ct C4 (|C c | 2 , |^-a;) 1 g T ®a;9| 7oCc / (|C c |a;) da;. 

Proof. The relative velocity and peculiar velocities after collision are expanded in 
a power series of e. For i € H, we have 

\C e -eCi\ = \C C \ - ejgj-a + e 2 6a + 0(e 3 ), 

^ = C l+ e^-a 4 -e 2 ^-a 5 + 0(e 3 ), 
C' c = \C c \u> + ea 5 + e 2 (±a 4 + a 6 ) + 0{e 3 ), 



with 61 = ^[iai 2 -(|^-a) 

Hence, the distribution functions arc found to be 



, a 4 = C c -|C c |u>, a 5 = d-^-CiUj, a e = biU. 



e 2 ^-dcJi(Ciya 5 +0(e 6 ), igH, 



/ e (C e ') = /e(|C |w) +e9 Co / c (|C c |a;).a5 + - £ ^^ Ce / e (|C e |u;):a 5 ®a 5 

+ e 2 d c J c (|C e |«). f^a 4 + a 6 ) + 0(e 3 ). 



Combining these equations, we obtain eqs. (|3.12|) - (|3.14|) after some lengthy calcu- 
lation. □ 

Theorem 13.21 admits three corollaries. First, we define the entropy produced at 
order e° in collisions between electrons and i heavy particles 



T ™ = -j ^(/e,/i)(Ce) In 



(27r) 3 / 2 n° 



fe (Co) 



dC , i £ H, 



where Q° = (27T TO^ksT /hp) 3, 72 is the translational partition function of electrons. 
Symbol hp stands for Planck's constant. The zero-order operator describes the 
relaxation of the electron populations towards an isotropic distribution function in 
the mean heavy-particle frame. 

Corollary 3.4. The zero-order collision operator Soiife > ft), i£H, vanishes when 
f c is an isotropic function of the velocity C e . Moreover, the zero-order entropy is 
nonnegative, that writes > 0, i £ H, and the inequality is an equality if and 
only if f e is an isotropic function of the velocity C c . 

Proof. If f e is an isotropic function of C e , we have f e (|C c |o;) = f c (C e ) for any 
<jj in the unit sphere, so that expression (|3.12[) implies that 2%{fe,fi) = 0. The 
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zero-order production entropy production reads 



^(iCd 2 ,^) 



/ e (|C e |u;)-/ e (C e ) 
(2n) 3 / 2 n° 



In 



Qe 



fe(C e ) 



d|C e |d^du;, 



and intertwining with u), 



c 

ICS, I 



|C C | 3 0(/ C (|C e |w), / c (C e ))d|Ce|d^du;. 



where Q(x,y) — (x — y)hx{x/y) is a nonnegative function. We then obtain that 
Tg if i 6 H, is nonnegative and equal to if and only if f e is isotropic in the mean 
heavy-particle frame. □ 

Corollary 3.5. The first-order collision operator 3ei(fe> fi)> i 6 H, vanishes when 
fi is an isotropic function of the velocity Ci . 

Proof. Expression (|3.13[) immediately yields that the integrand is an odd function of 
Ci, i € H, if f is isotropic in the mean heavy-particle frame, so that the first-order 
collision operator vanishes. □ 

Corollary 3.6. Considering f° = n c exp [-C e -C e / (2T e )] / (2irT c ) 3/2 and fP = 
rijm^ 2 exp [— mjCj-Ci/(2T/j)] / (27rT/j) 3 ^ 2 , i G H. the second-order collision oper- 
ator2 2 ei (f e °JP)(C e ), isH, reads 



(3-15) 3l(fe,f)(C c 



(T h -T c )^±f c °(C c )\C c 
mi J„ 



\c\ l 



)Q. 



a)/ 



Proof. A direct calculation of fi){C e ) given in (|3.14[) immediately yields 

expression (j3.15[) if f® and f® are the Maxwcll-Boltzmann distribution functions 
given in the assumptions of corollary 13. 61 □ 



3.4. Electron and heavy-particle collisional invariants. In the asymptotic 
limit e — > 0, the space of collisional invariants X defined in eq. (|2 .33[) splits into two 
subspaces naturally associated with our choice of scaling. 



(3.16) 



Definition 3.1. The space of scalar electron collisional invariants T c is spanned by 
the following elements 

| & = 1, 

Definition 3.2. The space of scalar heavy-particle collisional invariants Ih is 
spanned by the following elements 

4>h = 3 G H, 

^h +v = (miQ») ieW ve {1,2,3}, 

$h +4 = (l m «^"^)ieH> 

where symbol rP denotes the cardinality of the set H. 



(3.17) 
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The decoupling of the collision invariants is clearly identified in the proposed 
scaling. More precisely, the definition of the electron linearized collision operator 
(given in Section 3]) will involve the electron partial collision operator 3 ee and the 
mixed partial collision operators 3^, i G H, satisfying the following important 
property. 



Property 3.1. The partial collision operators S%i i S H, are orthogonal to the 
space of collisional invariants T c , i.e., ((y e ,3%)) e = for all I £ {1,2}. 

Proof. The projection of the collision operator 3^, i6H, onto tp l e , I € {1, 2} reads 
«^X}) e = ^ /^(|C c | 2 ,u;.^) |C c | 3 [/ (|C e |u;) - / e (C e )] # d|C e |du;d 

Intertwining a; with the projection {(^i,dci)) c is shown to vanish for all I 6 
{1,2}. □ 

We underline that the partial collision operators 3^, z £ H, are not orthogonal 
for the scalar product ((•,•))_ to the space spanned by the electron momentum. It 
is the reason why the vector C e does not belong to I c . In contrast, the definition of 
the heavy-particle linearized collision operator (given in Section 2]) will only involve 
the heavy-particle partial collision operators 3ij, hj £ H. 

Subsequently, using the newly defined collisional invariants, the orthogonality 
property 12.21 of the cross-collision operators can be rewritten 

(3.18) E^e.Zei» e = 0, h ,3he)) h = O, » € H, 

for mass conservation, 

(3.19) e]T« C ->^»c + i(^f +v M) h = 0, v £ {1,2,3}, 
for momentum conservation, and 

(3-20) «4 2 ,0ei» e + ((i>f + \2he)) h = 0, 

for energy conservation. This set of relations is essential since it corresponds to 
the conservation of mass, momentum, and energy in the electron heavy-particle 
interactions through the various orders in e of the Chapman-Enskog expansion. 

Then, the macroscopic properties are expressed as partial scalar products of the 
distribution functions and the new collisional invariants 

f = A. 

U/e,V>e»e = A^e, 



(3.21) 
and 



'{{kAi)) h = Pi, ieH, 

i-^i { U,i>f + l) h = 0, ^£{1,2,3}, 

AkAt +A )) h = Pne h . 

Symbol e c stands for the electron thermal energy per unit mass and eh, the heavy- 
particle thermal energy per unit mass. It is important to mention that the later 
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quantities are defined in the mean heavy-particle frame. Furthermore, the de- 
coupling of the collisional invariants is also consistent with the expression of the 
macroscopic properties. In particular, the electron momentum is not a collision in- 
variant in the proposed asymptotic limit, the electron mass flux is not constrained 
in the mean heavy-particle velocity referential. 

Moreover, translational temperatures are introduced as averaged thermal ener- 
gies in the mean heavy-particle frame. 

Definition 3.3. The electron and heavy-particle translational temperatures are 
given by 

(3-23) Tc = ^-((f c ,tt}) c , 

(3-24) Th = ^- h ({fh,j>f + %, 

where the heavy-particle number density reads nh — EjeH n r 

Consequently, the energy can be rewritten 

3 

((/e.^e))e = ^ n c T ^ 

for the electrons, and 

{{fh,4t +i )) h = \n h T h 

for the heavy particles. Moreover, it will be shown in Section 2] that both temper- 
atures are generally different. 

4. Chapman-Enskog method 

We resort to an Enskog expansion to derive an approximate solution to the 
Boltzmann equations (|2.29p - (|2.30|) by expanding the species distribution functions 
as 

(4.1) fc = / e °(l + e4k + e 2 4 2 + e 3 4 3 ) + ^(e 4 ), 

(4.2) fi =f?(l + e<k +£ 2 ^) + 0( e 3), zeH , 

and by imposing that the zero-order contributions and f® yield the local macro- 
scopic properties 

(4.3) ((f e °A = {{fe,tt)) e , le{l,2}, 

(4.4) ((J, , ft h )) h = ({f h , ft h )) h , I e {1, . . . , n H +4}. 

Hence, based upon the dimensional analysis of section |2~41 the electron Boltz- 
mann equation (I2.29[) appears to be 

(4.5) £ - 2 ^- 2 (/ c °) + e- 1 ®- 1 ^, 4 ) + ®°(fe°, 4 , €) + 4 , C €) 

= e- 2 r c 2 + e-H- 1 + t + s3l + 0(s 2 ), 
where the electron streaming operators read at successive orders 
^c~ 2 (/c°) = S bl q e (C c AB)-d c J c a , 
^\f^ck) = t 1 (/ c °)+9 e (^C c AB)-9c c (/ c %), 

@«\f?) = TF h C c -d x f^ + qe{^- h E' + 5 b0 C e AB)-d c J°, 
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$>°(f c °, <k , 4 2 ) = ^c(/c°, 4 ) + q e (5blC e ABydc e (fe°<t>e), 

^ e °(/ c °, 4 ) = d t f c ° + ^-C c -d x (f c % ) + v h -d x f e ° - {dcJ c %C c ):d x v h 
+ q e (S b0 M h v h AB + 5 b{ - 1) C c AB)-dcJ c ° 
+ qe{^E' + 6 b0 C e AB)-dc e (f e °<k), 

®l{f c \ 4 , 4 2 , = ®l (/e°, 4 , €) + « e («WaAB).0o.(/e°4 3 ). 

^0°, ^ , 4 2 ) = d t (f% ) + ^c c -a x (/ c °4 2 ) + v h -dM e % ) 

- M h ^-d c J e ° - (d Cc (f c °ck)®C c ):d x v h 

+ qe(6 b (-i)M h v h AB + 6 H _ 2) C e AB)-d c J° 

+ g e (5hoM ft u h AB + *6 ( _i)CUB).0 c .(/e O &) 

+ <7e(^£' + 6 b0 C e AB)-dc B {fe€ ), 

with the electric field expressed into the mean heavy-particle frame E' = E + 
S b iM?VhAB. The electron collision operators are given by 

for ease of readability in Sections 14. 2114. 7[ we strike through the collision operators 
that vanish when f® and i 6 H, are isotropic functions, 

37 1 = 3cc(/ c °4 , /e°) + 3ec(/ °, / c °4 ) 
J'EH 

t = M/cVc, / e °) + 3cc(/ c °4 , / c °4 ) + 3cc(/ c °, f e °<£) 

+ E 4(/c°4 2 , f?) + 3%l£dM)+ 3°, 

t = E 4-(/c ^ . ) + 3^%^+ 4(/ e °, ) + 4(/ e °, //), 

3c 1 = 3cc(/ c °4 3 , /c°) + 3cc(/ °4 2 , /c°4 ) + 3cc(/ °4 , /e°4 2 ) + 3cc(/ c °, fU!) 
+ E CC/eVe, J0>) + ^J^WT+ 31 

3c = E R'C/eVe, //4 ) + 3^(/ c °4 , + ^U^7jf)+ 4(/ e °4 , //fc ) 

j'eH . 

4 (/c°, J&°#) + 4 (/cV , //) + f e M , &4i ) + 2ei(/e°, //)}• 

Likewise, the heavy-particle Boltzmann equation (|2.30p is found to be 
(4.6) ^(f?)+e®l{f?,<k)=e-% 1 +f i +ed\+0{e 2 ), i e H, 
where the heavy-particle streaming operators read at successive orders 

= dtf" + {^-C l+ v h )-d x f? + + 6 bl C iA B).d c J° 

-M h ^-d c J? - {d c J l \C t ):d x v h , 
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®l(f°, <k ) = d t {f% ) + {^C l+ v h )-d x (fX ) 

+ %:6 b0 [(Q+Mhv^aB] -d c J t + + 5 bl C^B)-dcAf l %) 

The heavy-particle collision operators are given by 

37 1 = E i^fy+Mfi, 

t = e ^ ^ - //) + «° ) + iu^%rry + 3i 

d i =di(f,f e %) + 3l(f°,f e ) > 
j en 

3} = 3}M% , /cV ) + 3l(ff, /eVc 2 ) + dl(f% , / c °) + 3t(f, f e % ) + 

In the Chapman-Enskog method, the plasma is observed at successive orders of 
the e parameter equivalent to as many time scales. The micro- and macroscopic 
equations derived at each order are reviewed in Table [2] 

Table 2. Chapman-Enskog steps. 



Order Time Heavy particles 



Electrons 



>h 



Expression of 
Thermalization (T e 

Expression of f°, i £ H Equation for <fh 
Thermalization (Th) 

Equation for fa , i £ H 
Euler equations 



Zero-order momentum relation 
Equation for fa\ 

Zero-order drift-diffusion equations 
First-order momentum relation 



Navier-Stokes equations First-order drift-diffusion equations 



4.1. Order e~ 2 : electron thermalization. We resolve the electron Boltzmann 
equation (|4.5[) at order e~ 2 corresponding to the kinetic time scale t®. The electron 
population is shown to thermalize in the mean heavy-particle frame to a quasi- 
cquilibrium state described by a Maxwell-Boltzmann distribution function at tem- 
perature T c . In contrast, heavy particles do not exhibit any ensemble property at 
this order. 

Proposition 4.1. The zero-order electron distribution function f®, solution to 
eq. (|4.5[) at order s~ 2 , i.e., 2>~ 2 {f®) = 0~ 2 , that satisfies the scalar constraints 
(|4.3p is a Maxwell-Boltzmann distribution function at the electron temperature 



.3/2 



KINETIC THEORY OF PLASMAS: TRANSLATIONAL ENERGY 



23 



dC n = 0, 



Proof. Multiplying the equation @^ 2 (f°) = 3~ 2 by In [(27r) 3 / 2 n°/ c °/Q°] and inte- 
grating over dC c yields the entropy production 

T° + £ T e, +S biq J(C c ,BydaJ c °\n [(2*)»/Wjf /Q°" 

jGH J 

with T° cc = -/3co(/ °,/ °)(Co) In [(27r) 3 / 2 n°/ c °/Q0] dC . Using the equality 

a Cc / c °ln[(2^) 3 / 2 n / o 7^] = 9 Co {/ c °ln[(2^) 3 / 2 n /c7^] -/c°} and integrating 
by parts, the entropy production is found to be T° G + Yljen ^cj = 0. Moreover, a 
well-established derivation yields 



T° — — 



\ J n(/ B °j£ , / e °7ei)|Ce - C cl |a col du;dC dC cl > 0. 



Using coroilarv l3.41 we first obtain that > 0, i £ H, so that both terms T GC = 
and T Gi = 0, i G H. Then, corollary 13.41 implies that fP is isotropic in the mean 
heavy-particle frame. Seeing that T GO = 0, ln^, is thus a collisional invariant, i.e. 
is in the space T c . By using the macroscopic constraints, expression ()4.7j) is readily 
obtained. □ 

The choice of the referential velocity in which electrons thermalize will turns 
out to be crucial for the rest of the resolution. In the u = Vh frame, the quasi- 
equilibrium electron velocity distribution function is isotropic and the electrons 
follow the bulk movement associated with the hydrodynamic velocity of the mixture, 
following a physically relevant scenario. As already mentioned, the mean heavy- 
particle velocity does not depend on the small e parameter while still being 
close to the actual hydrodynamic velocity v of the entire mixture; this property is 
essential in order to conduct a rigourous multiscale analysis in the framework of the 
present Chapman-Enskog expansion. The relevance of such a choice of referential 
will be thoroughly investigated in section [ 



4.2. Order e" 1 : heavy-particle thermalization. We resolve the heavy-particle 
Boltzmann equation (|4.6|) at order e -1 corresponding to the kinetic time scale 
The heavy-particle population is shown to thermalize in the mean heavy-particle 
frame to a quasi-equilibrium state described by a Maxwell-Boltzmann distribution 
function at temperature Th- 

Proposition 4.2. Considering f® given by eq. (|4.7p . the zero- order family of heavy- 
particle distribution functions f£ solution to eq. (|4.6p at order e -1 , i.e., 3^ = 0, 
! 6 H, that satisfies the scalar constraints (|4.4j) is a family of Maxwell-Boltzmann 
distribution functions at the heavy-particle temperature 

Proof. As the zero-order electron distribution function f? is isotropic in the mean 
heavy-particle frame, corollary 13.11 yields that the heavy-particle Boltzmann equa- 
tion (|4.6[) reads at order e _1 

jen 

After some classical algebra [2] , we obtain expression (|4.8[) for the zero-order heavy- 
particle distribution functions. □ 
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Thus, propositions 14.11 and 14.21 allow for electron and heavy particles quasi- 
cquilibrium states to be obtained at different temperatures. 

4.3. Order e _1 : electron momentum relation. We conduct the resolution and 
derive a momentum relation based on the electron Boltzmann equation (14. 5p at 
order e _1 corresponding to the kinetic time scale t^. We then emphasize an original 
property of the Chapman- Enskog expansion at this order associated with both the 
absence of momentum constraint in eq. 13.211 and the multiscale analysis. 

With the previously obtained Maxwell-Boltzmann electron distribution function, 
we first define the electron linearized collision operator in the case. 

Definition 4.1. The electron linearized collision operator 5F C reads 

?e(4 ) = ~ [3ce(/ c °4 , / c °) + 3 cc (f c ° , f c % ) + £ <k , //) 



where f° is given by eq. (|4.T[) and f° by eq. 

The kernel of this operator is given in the following property. 

Property 4.1. The kernel of the linearized collision operator J c is the space of 
scalar electron collisional invariants T e . 

Proof. The linearized collision operator 3" e is rewritten in the form 

%(<k ) = - J /el (<& + <4'l - 4 - 4l ) |C e - Coi keel du;dC cl 

-J2 n i I ^(\C c \ 2 ,u-^\C c \(ck(\C c \u:)-ck(C c ))du:. 

We then obtain that the space I e is in the kernel of 3e- Reciprocally, if 9^(4 ) = 0, 
multiplying 3^ e (4k ) by f®4k and integrating over dC e yields 

\ j /c7el (<& + <4'l ~ <k ~ iOVe - Celkeel du;dC c dC cl 



2 ^ 



j y do, (|C | 2 , U;-|g|) |C C |/ C °(4 (|Ce|u;) - 4 (C c )) 2 du;dCe = 0, 



so that 4k is in the space X e . □ 

Based on corollaries 13.41 and 13.51 the electron Boltzmann equation (|4.5[) is found 
to be at order e -1 

(4.9) f c °?M + S bl q e dc B (fe<k)-C e AB = -@-\f c °), 

with the constraints 

(4-10) ({f e %,i>l)) e = 0, I €{1,2}. 

The terms dc (fo4h )-C c /\B and ^~ 1 (f°) are orthogonal to the kernel of 5" e for 
the scalar product ((■, •)) . Consequently, no macroscopic conservation equations of 
mass and energy are derived at this order. 

In fact, for any value of to, defining the shifted Maxwell-Boltzmann distribution 

(4-11) r° = n c (^r) 3/ exp (~(G, - eM h wf 
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we can expand it as a function of e 

)+0(e 3 ), 

which still yields, at leading order, the same distribution as defined in cq. (I4.7|) . 
We then realize that the Chapman-Enskog expansion can be rewritten in a different 
way at this order : 

(4.13) / e °(l + e4 + £ 2 4 2 ) = r°(l + + + 0(s 3 ), 

with 
(4.14) 

-L c 

It is interesting to notice that, whatever the choice of w, the part of the hydro- 
dynamic velocity of the full mixture 

(p h + e 2 p c )M h v = p h M h v h + e 2 p c v c , 

associated with the electrons p e v e will be splitted into two parts at the same order 
of the multiscale expansion 

v c = M h v h + V C + 0(e) = M h (v h + w) + V c w + 0(e), V c w = V c - M h w, 

with p c V c w = J C c f™®4>™&C c . Thus, as opposed to the classical expansion, since 
no momentum constraint is to be found for the electrons, the definition of the mix- 
ture hydrodynamic velocity does not allow to uniquely define the electron diffusion 
velocities. In any case, the hydrodynamic velocity of the mixture is Vh at order e . 
It is then necessary to properly delineate the possible choices for the w velocity, 
which should not be confused with a change of referential, since it only influences 
the electron Chapman-Enskog expansion. 

Lemma 4.1. In the chosen frame of reference, any velocity w leads to a new 
definition of <t>™ for which property ^. 1\ is preserved and thus leads to an equivalent 
resolubility condition for 4>™ as for 4 . Moreover, the resolution of 4>™ is completely 
equivalent to the resolution of 4 . 

Proof. It is sufficient to note that the difference 6$f = </>™ — 4k = —MhC c -w/T c is 
orthogonal to the collisional invariants {(fcd&^iWc^c = 0, Z G {1, 2}. □ 

It is interesting to note that for our choice of moving frame u — Vh, the electron 
thermalization naturally occurs in the "appropriate" referential in close connection 
to the physics of the problem and there is no need to use the abovementioned 
property in order to conduct the resolution at order e _1 , therefore, we take w = 
in the following. We will also have to check the validity of such a strategy at higher 
orders; we will come back to this point in Section [4.51 

As mentioned earlier, the partial collision operators deii i£H, are not orthogonal 
to the space spanned by the vector C c . However, an electron momentum relation 
is obtained by projecting eq. (|4.9p onto this space. First, the electron pressure, 
diffusion velocity, mean velocity, conduction current density in the mean heavy- 
particle velocity frame, and conduction current density in the inertial frame are 



(4.12) 



/°( 1 + ef^ + e*§ 



(c c - w y 



iw2 



M_h 



-(C c -w) 



Ml 

2T n 



—ww 



(C c -w) 2 
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defined as 

(4.15) p c = n c T c , 

(4.16) V e = — [ C e f e % dC c , v c = M h v h + V c , 

n e J 

(4.17) J c = n c q e V e , j c = n c q e v c . 
Then, we have the following proposition. 

Proposition 4.3. Considering / G ° given by eq. (|4.7[) and f° , i E H, by eq. (|4 

the zero-order momentum exchanged between electrons and heavy particles reads 

(4.18) J2 «3c,(/c°4 , //), C* )> c = -Lfl^pe - -rr-E - 5 b ij c *B. 

jGH h h 

Proof. Equation (|4.9p is projected onto the space spanned by the vector C c 

-{(f e °%(ck),C e )) e = 0-\C c )) c + 5 bl q e (d Ca {f^)-{C c AB),C c )) c . 

Then, eq. (14. 18)) is readily established by simplifying the left-hand-side by means 
of eq. (|2.36p . ((C c , 3co)) = 0, at order e and by integrating by parts the right-hand- 
side. □ 

The zero-order momentum exchanged between electrons and heavy particles is 
thus expressed in terms of the electron pressure and electric force. In addition, the 
following lemma allows for the momentum exchanged between heavy particles and 
electrons to be calculated at order zero. 

Lemma 4.2. Considering f® given by eq. (|4.7[) and f®, i G H, by eq. (|4.8|l . the net 

zero-order momentum exchanged between electrons and heavy particles vanishes, 
i.e. 

(4-19) ({31 U°, f c % ) , ^ h+v )) h + ^({3% (f c % , , c eu )) e = 0, 
for {1,2, 3}. J£R 

Proof. Equation (|4.19[) is derived from eq. (|3.19[) at order e 2 based on corollaries 
OGL3EH andE31 □ 

Moreover, the zero-order momentum exchanged between heavy particles and 
electrons can be directly calculated after introducing the average force of an electron 
acting on a heavy particle i given by 

(4-20) F ie = J Q^dTel 2 ) I7cl7 c Zc°(7c)<M7c) d le , i E H. 

Lemma 4.3. Considering f® given by eq. (|4. T[) and f®, i € H, by eq. (|4. 8[) . the 

zero-order momentum exchanged between heavy particles and electrons reads 

(4.2i) ((si a°, f c °<k ) , r h h+v )) h = £ n ^ 

for ve {1,2,3}. 

We will see that the average forces F ie , i £ H, contribute to the heavy-particle 
diffusion driving forces and, in particular, yielding to anisotropic diffusion velocities 
for the heavy particles in the case 6=1. 
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4.4. Order e°: heavy-particle Euler equations. We derive Euler equations 
based on the heavy-particle Boltzmann equation (|4.6j) at order e° corresponding to 
the macroscopic time scale t . First, a linearized collision operator is introduced 
for heavy-particles. 



Definition 4.2. The linearized collision operator 3^ = (3*i)i e H reads 

i G H, 



where f®, i G H, is given by eq. (|4.8[) . for a family = (0i)i S H- 

The first non-vanishing term of the partial collision operator 3 he is not included 
in the linearized collision operator since it does not exhibit any property of orthogo- 
nality to lh for the scalar product ((■,■)) h - The kernel of 3^ is given in the following 
property, the proof of which is omitted since it is a well-established result [25] . 

Property 4.2. The kernel of the linearized collision operator 3tj is the space of 
scalar collisional invariants Th- 

Furthermore, we define the heavy-particle pressure, ph = n^Xh, the mixture 
pressure, p = p c + Ph, the heavy-particle charge, nhqh = J2j£H n j1j^ the mixture 
charge, nq = n c q e + nhqh, and the total current density Iq = UhQh v h + n c q e v e /Mh- 
The energy exchanged between heavy particles and electrons reads at order zero 

(4-22) AE° = ({3UgJ c °)^: h+i )) h - 

This quantity is of the order of the thermal energy divided by the macroscopic time 
scale, n°kBT°/t . A more accurate expression is calculated by means of corollary 



(4.23) A£g = (T e -r fc )£n i ^ 



11 
m,- 



Then, the heavy-particle Euler equations are derived in the following proposition. 
Proposition 4.4. If <f>h is a solution to eq. (14. 6|) at orders , i.e. 
(4-24) J?Wh) = -2?(l?)+3l iGH, 

where f® is given by eq. gZ3), //>, i£H, by eq. and 4 by eqs. (|i^ |) -(|4TTD |) . 

and if f^4>h = {fi<k )iGH satisfies the constraints 

(4.25) «/h%,&» fc = 0, le{l,...,nP+4}, 

i/ien, the zero-order conservation equations of heavy-particle mass, momentum, and 
energy read 

(4.26) d tPi + d x -{ Pi v h )=Q, i G H, 

(4.27) d t (p h v h ) + d x -{p h v h ®v h + jjppT) = jpnqE + S bl I AB, 

h h 

(4.28) 9t(/3 h e h ) + d x - (p h e h v h ) = - Ph d x -v h + AE°. 

Proof. Fredholm's alternative [30] represents the solvability condition of eq. (|4.24|) 
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/ G {1, . . . , rfi + 4}. Integrating by parts the left-hand-side and simplifying the 
right-hand-side based on theorem 13.11 and corollary 13. 2\ one obtains eqs. (|4.26l) . 
(|4.28p , and the following momentum conservation equation 
(4.29) 

- M hPh ^ - ^d xPh + j^n h q h E' + ({31 (4°,/ e °4) , (r h h+ l,e{i,2, 3} )) h = 0. 

Simplifying the latter equation by means of the heavy-particle mass conservation 
equation dtPh + 9 x -(ph,Vh) — and lemma [4T21 yields eq. (|4.27[) . □ 

4.5. Order e°: zero-order electron drift-diffusion equations. We derive zero- 
order electron drift-diffusion equations and a momentum relation based on the 
electron Boltzmann equation (|4.5[) at order e° corresponding to the macroscopic 
time scale t°. We also prove, at this order of the resolution, that any nonzero shift 
introduced at the previous order leads to a series of difficulties at the present order. 
It thus demonstrates that the initial choice of referential leads to a quite natural 
resolution at sucessive orders. 

With the previously obtained Maxwell-Boltzmann electron distribution function 
in eq. (|4.7p we introduce the electron heat flux 

(4.30) q e = [ lc c -C c CJ c °<kdC c . 



2 

The energy exchanged between electrons and heavy particles reads at order zero 

(4.31) A£; o = ^((4(/ c ,^),v3 o 2 )) c . 

JGH 

The latter expression is calculated by means of eq. (|3.20| at order e 2 

(4.32) AE° + AE° = 0, 

where AE° is given by eq. (|4.23|) . Then, the zero-order electron drift-diffusion 
equations are derived in the following proposition. 

Proposition 4.5. If 4>% is a solution to eq. (|4.5|) at orders , i.e. 

(4.33) f?3M) + s bl q e d Co (/ c v c 2 )-c ab = -# e °(/ e °, 4 ) + a C c(/ c °4 , /e°4 ) + 21 

where f° is given by eq. (fOTl . f° , i G H, by eq. (|4.8[l . (jh by eqs. I|4.9 |l - (|4.10p . and 
4>i , i G H by eqs. (|4.24| - ()4.25p . and if f?(/> 2 satisfies the constraints 

(4-34) ((/ o V e 2 ,^)) o = 0, le{l,2}, 

then, the zero-order conservation equations of electron mass and energy read 

(4.35) d tPc + d x -(p c v h + ^p c V c ) = 0, 

(4.36) d t { Pc e c ) + 8 X - { Pc e c v h ) = -p c d x -v h - jr h d x -q c + j^J c -E' + AE°. 

Proof. Fredholm's alternative [3U] represents the solvability condition of eq. (|4.33p 

0li>i)) e = {§lAi)) e , ^{1,2}. 

Integrating by parts the left-hand-side and simplifying the right-hand-side based 
on theorem 13.21 and corollary 13.51 yields eqs. (|4.35|) and (|4.36[) . □ 
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Lemma 4.4. In the chosen frame of reference, any velocity w leads to a new defi- 
nition of <j)™ 2 in eq. (I4.14| . for which property ^. 1\ is preserved, and thus leads to an 
equivalent resolubility condition for <f>™ 2 as for </> 2 . However, the resolution of 4>™ 2 
is not equivalent to the resolution of 4 2 ■' in particular, the expansion corresponding 
to w =/: yields a non standard Chapman- Enskog expansion where the second- order 
distribution perturbation does not satisfy the scalar constraints (|4.34|) . 

Proof. The difference between (j)™ 2 and (f> 2 reads (|4. 14[) 

= c 2 - € - -^(a-w)ck + H [„.«, + 

The projection of 6(f)™ 2 onto the collisional invariants is given by 

((/c^C 2 ,^ = ^n c w(M h w - V e ), 
((/ C (5C 2 > C 2 )) C = M h w(2M h n c w - lq e ). 

The difference b~4>™ 2 is then orthogonal to the collisional invariants if and only if 
w = 0. To conclude, the resolution of </>™ 2 yields a linearized Boltzmann equation 
where the right-member is orthogonal to the collisional invariants — a direct calcu- 
lation shows that 3 (S>t>™ 2 ) +5biq e dc (S4>™ 2 )-C aB is orthogonal to the collisional 
invariants — whereas the scalar constraints on the unknown function <f>™ 2 are not 
zero. □ 

Consequently for the reasons invoked so far, we will not try to shift the center 
of the Maxwell-Boltzmann distribution for electrons and stick with w = at any 
order. 

Then, we define the electron viscous tensor, second-order electron diffusion ve- 
locity, and second-order current density 

(4.37) H ° " / Cc ® Cc ^ dCc ' 

(4.38) V 2 = — [ C c / C °4 2 dC c , 

n c J 

(4.39) J 2 = n c q e V 2 . 

A first-order electron momentum relation is given in the following proposition. 

Proposition 4.6. Considering f® given by eq. (|4.7|l . f°, i G H, by eq. (|4.8|) . 4k by 
eqs. flUSJ-flHBJ), 4\ , i G H, by eqs. (|433j) - ([Qgjl . and 4 2 by eqs. (j4T3^]) - (|4~34l) . the 

first-order momentum exchanged between electrons and heavy particles reads 

(4.40) Y, ((4-U°4 2 , //), Cc}) + «3 C °, C » = ^d x -U c - (5 bo j c + S bl J 2 ) AB. 
Proof. Equation (|4.33[) is projected onto the space spanned by the vector C c 

-((/e°?e(4 2 ). C e » e + ({1°, C e )) e = 9°(/ e °, & ) + hide ({d Ce {£<£)-C kB, C e » e . 

Then, eq. (|4.40p is readily established by simplifying the left-hand-side by means of 
eq. (|2.36|) . ((C e , 3ee)) e — 0, at order e 2 and by integrating by parts the right-hand- 
side. □ 
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The first-order momentum exchanged between electrons and heavy particles is 
thus expressed in terms of the electron viscous tensor and electric force. Besides, 
the following lemma allows for the momentum exchanged between heavy particles 
and electrons to be calculated at order s. 

Lemma 4.5. Considering f® given by eq. (|4.7[) . fP, i G H, by eq. (|4.8|) . 4k by 
eqs. (|4T9 ]) -(|430 |) , ^,jgH, by eqs. (|4TM|) - (|4^5|) . and $ by eqs. (|433>(l434j) . the 

net first- order momentum exchanged between electrons and heavy particles vanishes, 
i.e., 

(4.41) ((3l^ h+ n) h + E^ei(/e°4 2 , C eu )) e + ((2°, C e „» e = 0, v e {1, 2, 3}. 

Proof. Equation (|4.41 1) is derived from eq. (|3 . 19[) at order e 2 based on corollaries 
1013^61 □ 

4.6. Order e: heavy-particle Navier-Stokes equations. We derive Navier- 
Stokes equations based on the heavy-particle Boltzmann equation (|4.6[) at order e. 
First, we introduce the diffusion velocity and mean velocity of species ieH, 

(4.42) Vi = - [ CifP&dd, v i = v h + -faV i , i6H, 
the heavy-particle viscous tensor, 

(4.43) n, = ]T / n,,C, CJ,'o,<\C r 

the second-order electron mean velocity, 

(4.44) v 2 c = M h v h + V c +eV c 2 , 
the heavy-particle heat flux, 

(4.45) q h = [ ^"'^yCXhrOi dC* , 

the heavy-particle conduction current density in the mean heavy-particle veloc- 
ity frame, the heavy-particle conduction current density in the inertial frame, the 
second-order electron conduction current density in the inertial frame, and the total 
current density, 

(4.46) J h = ^n q V v j h = ^n q ] v j , £ = n c q e v 2 c , I = j h + ^j-j 2 . 

Furthermore, we define the energy exchanged between heavy particles and electrons 
reads at order e 

(4.4?) aei = {(di(fh<k , f c %),r h h+4 )) h + (tfUfhi 

+ ((3UfH<kj c °),r h h+i )) h - 



The first term can be calculated by means of theorem 13.11 

(4.48) ((aL(/, 4,/ c 4)^r +4 )), = Y, n i v i' F i« 

and the two other terms will vanish. Then, we establish the following lemma used 
in the derivation of the heavy-particle Navier-Stokes equations. 
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Lemma 4.6. Considering f° given by eq. (|4.7|) . f®, i G H, by eq. (|4.8jl . fa by 
eqs. (f4T9 ]) - ff4~T0|) . fa, i G H, by eqs. (f4~^4|) - (j4~25|) . and 4 2 6y eqs. <f4^3^ - f|4734|) . i/ie 
mass exchanged at order e between heavy particles and electrons vanishes, i.e., 

(4.49) = 0, lG{l,...,n h }. 

Proof. Equation (|4.49[) is readily derived from eq. (I3.18| at order e 3 . □ 
Proposition 4.7. If ' fa\ is a solution to eq. (|2 . 30[) at order e 1 , i.e. 

(4.50) ffli (fa 2 , ) = -@\ (f, fa) + Y, 3ij (f°<k , $<k )+3h »£H, 

where / e ° is given by eq. (|4~7|) . //\ i £ H, by eq. (|4T5|) . &?/ egs. (|43 1) - (j4~TU)) . ^ , 
j £ H, 61/ egs. g^l-g^U), and # e<?s. ^M~^M, and if = {tffa 2 ^ 
satisfies the constraints 

(4.51) ((.f Q h tiJ l h )) h =0, ;e{l,...,n H + 4}, 

then, the first- order conservation equations of heavy-particle mass, momentum, and 
energy read 

(4.52) d tPi + d x -{ Pi v h + -^piVi ) = 0, iGH, 

(4.53) dt(phVh) + d a .-(p h v h ®Vh + aTJP 1 ) = -jjfjfliB-Cnfc + n c ) + jj?nqE 

+ [S b0 I + 6 bl l]AB, 

(4.54) dt(pheh) + d x -(p h e h v h ) = -(p/J + eU h ):d x v h - j^-d x -q h + j^J h -E' 

+ AE° h +sAEl 

Proof. The Chapman-Enskog method allows for the following conservation equa- 
tions to be derived 

& h )) h + e((9lft h )) h = ((3li> l h )) h + e((3l^ h )) h , 

I G {1, . . . , n H + 4}. Integrating by parts the left-hand-side and simplifying the 
right-hand-side based on the proof of heavy-particle Euler eqs. (|4.26[) - (|4.28|) , propo- 
sition rj3J lemma 1431 and lemma 1431 one obtains eqs. (|4.52|) - (|4.54|) . 

□ 

Remark 4.1. When one single type of heavy particles is considered, the first-order 
energy exchange term, heavy-particle diffusion velocities, and conduction currrent 
degenerate, Ai?^ = 0, V i = 0, i G H, Jh = 0, the total current is simplified 
as well, I = nqvh + n c q e V c /Mh. Therefore, we retrieve the formalism of Degond 
and Lucquin. In such a case, the Navier-Stokes system can be coupled to the 
system of drift-diffusion equations for the electrons obtained at order e° in the 
previous section. Since no energy exchange occurs at order e , there is no need 
to resolve the electrons at order e 1 to obtain a conservative model which insures 
the positivity of the entropy production. However, this oversimplified case hides 
the details of the complex interaction between the electron and the heavy-particle 
mixture which is exhibited by the previous system of conservation eqs. (|4.52l) - (|4.54l) . 
For a multicomponent mixture of heavy particles, thus, we have to extend the model 
obtained so far for the electron one order further, as done in the following section. 
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4.7. Order e: first-order electron drift-diffusion equations . We derive first- 
order electron drift-diffusion equations based on the electron Boltzmann equa- 
tion (|4.5[) at order e 1 . 

Then, we define the second-order electron heat flux 

(4-55) q 2 c = J iCe-CeCe/ > c 2 dC c . 

The energy exchanged between electrons and heavy particles at order e is caculated 
by means of eq. (13.20|) at order e 3 

(4.56) AE X C + AEl = 0, 

where AE^ is given by eq. (|4.47[) . Moreover, we establish the following lemma used 
in the derivation of the first-order electron drift-diffusion equations. 

Lemma 4.7. Considering f® given by eq. (|4.7[) . f®, i G H, by eq. (|4.8|) . 4 by 
eqs. (|4T9 ]I - (|4TTQ|) . fa , i g H, by eqs. (|i3gj) . and $ by eqs. (|4l%j) - (j434) . the 

mass exchanged at order e between electrons and heavy particles vanishes, i.e., 

(4-57) = 0- 

Proof. Equation (|4.57[) is readily derived from eq. (|3 . 1 8[) at order e 3 . □ 
Proposition 4.8. If faf is a solution to eq. (|4.5p at order e 1 , i.e. 

(4.58) / c °? c (4 3 ) + 8 bl q e d c . (/c°4 3 ) C c f\B = -®\ (f c °, fa,^) + M/eVe, / c °4 ) 

+ aec(/ c 4,/c°4 2 )+^, 

where f® is given by eq. (|47f|l . fP, i e R, by eq. JO), 4 fey eqs. (|43D- (|4~TU1) . 4, 
i g FL 6?/ egs. (14^24) - (14^5)) . 4 2 by eqs. g^-g^l, and 4?, i <EH, by eqs. (j43fj|) - 
(|4.5fp . and if f?(f> 3 satisfies the constraints 

(4-59) ((/ o >3,^)) e = 0, Jg{l,2}, 

i/ien 7 i/ie first-order conservation equations of electron mass and energy read 



(4.60) d tPc + d x - Pc (v h + j±-(V e + eV e 2 )) 



= 0, 



(4.6f ) dt(p c e c ) + 9 X - (p e e e v h ) = -p c d x -v h - j^d x - (q c + eq 2 c ) 

+ j^{J c + eJ 2 ) -E' + 5 bo eM h J e -v h AB + A£ c ° + eAE l e . 

Proof. The Chapman-Enskog method allows for the following conservation equa- 
tions to be derived 

+ = + e((3l&)) e , I G {1, 2}. 

Integrating by parts the left-hand-side and simplifying the right-hand-side based 
on lemma l4~7l one obtains eqs. (|4.60|) - (j4.61|) . The added terms associated with 
the perturbation of the mean electron velocity at order e 1 do not bring in any 
contribution to the conservation equations. □ 

Before switching to Section [5] in order to evaluate the expression of the transport 
fluxes, we briefly come back to question of the influence of the choice of the initial 
frame of reference. 
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4.8. About the necessity of conducting the expansion in the Vh frame. 

As mentioned earlier, the mean heavy-particle velocity frame is not commonly 
adopted in the literature to conduct the Chapman-Enskog expansion. We have 
already underlined that the natural choice of the hydrodynamic velocity frame 
is not appropriate insofar as the global hydrodynamic velocity v depends on the 
parameter e. Besides, the choice of the inertial frame provides a vanishing mean 
velocity of the electrons. In fact Degond and Lucquin [20] and Lucquin [361 EI] 
reach such a conclusion. However, since the expansion of the collision operators in 
terms of e depends on the initial choice of referential (see remark l3~Tj) and since the 
choice of the inertial frame prevents some terms from vanishing (such as ■ (fci fj*))? 
we will first show that these authors compensate the presence of non-zero terms 
in the integro-differential equations by the help of the w velocity on the electron 
introduced in Section 14.31 This is affordable for the resolution oi 4k, as proved 
in the following. We then investigate if such an approach can be extended to the 
resolution of 4e- 

Let us review the Chapman-Enskog expansion in a general frame. Considering 
a frame moving with the velocity it, the peculiar velocities are given by 

(4.62) C™ = c c sM h u, C? = Ci- M h u, i e H. 

The space of scalar electron collisional invariants T£ is spanned by the following 
elements 

fa 1 = 1, 

^ = \c:-c:, 

the space of scalar heavy-particle collisional invariants Z" by 

{C J = (mi%). eH , j e H, 

ft'**" = M?) ieH , ^G {1,2,3}, 

c« H+4 = (|m iC r-q% H , 

and the macroscopic properties are expressed as partial scalar products of the dis- 
tribution functions and the collisional invariants 

((f a U ,W'% - P., 

and 

W&X = Pi. *GH, 

((fh,i>h nH+U )) h = PhM h (v hv - u v ), v e {1, 2, 3}, 

UT^H nU+4 )) h = ln h T h + \Ml Ph {u-v h ).{u-v h ). 

Similarly to the low Mach number approximation for neutral gases, we decouple for 
the electrons the thermal energy from the mixture kinetic energy. 

Then, we reword the Chapman-Enskog expansion of Section [4] in the u frame. 
First, let us now formulate two propositions from the begining of this section. 

Proposition 4.9 (Order e~ 2 : electron thermalization). The zero-order electron 
distribution function f™° , solution to eq. (|4.5p at order e~ 2 , i.e., ^ c ~ 2 (/ c 1 ' ) = 
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2~ 2 , that satisfies the scalar constraints (if™ , W c )) c = {(f™,V e )) c , I £ {1, 2} ^ is a 
Maxwell- Boltzmann distribution junction at the electron temperature 

1 \ 3/2 / 1 



(4.63) = exp^-Cr-^ 

Proof. The proof is identical to the one of proposition 14.11 □ 

Proposition 4.10 (Order £ _1 : heavy particle thermalization). Considering f™° 
given by eq. (|4.63|) . the zero-order family of heavy-particle distribution functions 
f™° solution to eq. (|4.6p at order e~ x , i.e., = 0, i G H, that satisfies the scalar 

constraints {{fh°,w h )) h = ((/"> ' e {!> • ■ ■ > nIi +4}, is a family of Maxwell- 

Boltzmann distribution functions at the heavy-particle temperature 

(4.64) 

3/2 



(S) ™p(-^[ C * U - M hiv h -u)]-[C?-M h (v h -u)}), zgH. 



y0 _ 

ij 



Proof. Multiplying the equation ^r 1 = by In [(2Tr) 3 / 2 n° f^°(m^Q° h )], where the 
heavy-particle translational partition function reads Q? = (27rm°kBT°/hp) 3 / 2 , in- 
tegrating over dCi, and summing over the heavy particles yields the entropy pro- 
duction 

i,j£H i6H 

where the partial entropy production terms read 

J a^Q - Ci | n(r°7/°', r°//°) dQdqdo;, i, j G H, 
T° = ^ / d c jr°\n da- J Q\l\hf) |7el7 c f c u ° d 7c , < £ H. 

Integrating by parts, the terms T° e , j £ H, vanish and we conclude as usual. □ 

At this step, two properties appear: the electron thermalization takes place in 
any velocity frame, whereas the zero-order heavy particle distribution functions do 
not depend on the selected frame. Indeed, we clearly have f u0 = ff, i G H, for all 
velocity it. 

Considering then the Boltzmann equation at order e _1 , the electron first per- 
turbation (j)™ satisfies the linearized Boltzmann equation 

(4.65) ? c (4") + 5 bl q e d Cc A€>C^B = --L#-i(/«0) 

Jc 

+T. n ^^h\c:?)\c:\{v h -u)-c:, 

ieH e 

with the constraints 

(4-66) ((fCf'S^o, l£{l,2}. 

The right-member of eq. (|4.65[) satisfies the constraint to be orthogonal to the 
collisional invariants, that is the solvability condition. Moreover, in order to avoid 
treating the newly introduced term in the integro-differential equation, one can use 
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the absence of momentum constraints on the electron distribution function and 
introduce a velocity shift w = Vh — u and notice that [36] 

(4.67) %{c:) = -Y J n l Q^{\c:\ 2 )\c:\c:, 

we thus obtain that the conduction of the Chapman-Enskog expansion in the u 
frame is equivalent to that in the frame with 

± c 

As already mentioned in Section 14.31 the electron velocity v c can be splitted into 
two parts at the same order of the multiscale expansion v c = M^u + V c u + 0(e), 
with V c u = V c + Mhivh — u). We have thus provided a nice interpretation of the 
algebra proposed in Lucquin [36j where the use of w = Vh allows to eliminate the 
presence of the term X)jeH 3lj(f^°, fj 10 ) hi the integro-differential equation for 0™ 
obtained when working in the inertial frame u = 0. 

As a conclusion, it amounts to "coming back" into the mean heavy-particle veloc- 
ity referential. Let us underline at this point, that the set of equations obtained for 
the heavy-particle Euler equations coupled to the zero-order electron drift-diffusion 
equations is identical to the set obtained in Lucquin [36 . While still equivalent 
at this order of the expansion to our study and yielding the same macroscopic 
equations, it leads to an artificial complexity. This is a first step in the justifica- 
tion of the choices made in Section 12.41 in terms of the referential and associated 
simplified algebra. However, at order e, which yields heavy-particle Navier-Stokes 
equations coupled to first-order electron drift-diffusion equations, we realize that 
such a compensation used through the velocity shift w has a nasty influence on the 
structure of the expansion at the next order (see Lemma 14. 4j) and hence will make 
the resolution of </> c 2 difficult. Concerning the heavy-particle Boltzmann equation 
at order e°, the first-order perturbation functions </>", i £ H, satisfy too eq. (|4.24[) . 
that implies that 0™ = <fo , i £ H. 

Finally, the second-order electron perturbation 0™ 2 satisfies 

(4.68) J C (C 2 ) + 5 biqe d Ca {€ 2 YC c *B = * 4«) 

Jc 

Jc Jc J£H 

with the constraints 

(4-69) <</c°C 2 ,C''>} c = 0, ZS{1,2}. 

In the general case, the right-member of eq. (|4.68[) is not orthogonal to the collisional 
invariants, and the solvability condition is not satisfied. More precisely, the two first 
terms J>°(/ e u0 , O// c u0 and dcc(f c u0 ^, /"V")//" satisfy the constraints for all 
velocity u, whereas a lengthly calculation yields that the projection of the last term 
over the collisional invariants if}™' 1 , I £ {1, 2}, does not vanish, except for u = Vh- 

As a conclusion on this matter, if the velocity u is not equal to v^, that is to say 
if the selected frame is not the mean-heavy particle velocity frame, the structure 
of the expansion of the collision operators in terms of e prevents from solving for 
<4 m2 . Moreover, the "trick" used in Degond and Lucquin at previous order for 
4k also brings in additional difficulties at subsequent order through the residual 
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terms associated with w. The choice of the mean heavy-particle velocity frame not 
only yields a simplified algebra — the parity properties of the distribution functions 
are in accordance to the expression of the multiscale collisional operators — but 
is mandatory in order to obtain the well-posed second-order problem for electron 
distribution, second-order electron diffusion velocity, and heat flux being crucial to 
balance the global energy equation. This statement fully justifies the choices made 
in this contribution. 

5. Transport coefficients 

In this section, we investigate the electron and heavy-particle perturbation func- 
tions in order to obtain the expressions of the transport fluxes. We only treat the 
case 6=1 corresponding to a strong magnetization (inducing anisotropic trans- 
port coefficients), the case b = corresponding to a weak magnetization will be 
investigated in a forthcoming publication. 

5.1. Extra notations for anisotropy. We introduce some extra notations in 
order to conveniently express the solution to the Boltzmann equation in the presence 
of a strong magnetic field. First, we define a unitary vector for the magnetic field 
23 = B/\B\ and also three direction matrices 

M^SaS, M x =II-£®23, M = 

so that we have for any vector x in three dimensions 

cc" = M"cc = cc-23 23, x 1 - = M ± x = x - x "B 23, cc° = M e x = Sax. 

In the {x,"B) plane, the vector cc" is the component of x parallel to the magnetic 
field and x 1 - its component perpendicular to the magnetic field. Thus, we have 
x = cc" + x ± . The vector cc lies in the direction transverse to the (cc,23) plane. 
The three vectors cc", cc x , and cc are then mutually orthogonal. We will show 
that the transport properties are anisotropic. In the weak magnetization limit, the 
transport properties are identical in the parallel and perpendicular directions and 
vanish in the transverse direction. 

5.2. First-order electron perturbation function. The first-order perturbation 
function 4k is a solution to eq. (|4.9p 

(5.1) 'J c {(k)+qedcS ( k)-C^B = * c , 

and satisfies the constraints (|4.10[) . where ^f e is given by the expression ^> c — 
~^e 1 (fe)/fe aim fe by eq. (|4.7|) . After some algebra based on the expression of 
f®, the quantity *f! e is transformed into 

(5.2) * e = -p e *^-d e -*^-S a / ' 




where the electron diffusion driving force d c is defined by the relation 

(5.3) d c - -d xPe - ^E', 

Pc Pc 

and with 

(5.4) = — !— C c , *^ = — (~T e - -C c -C c )c e . 
v ' c M h p e ' M h \2 2 J 
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The right-hand-side of eq. (15. ip does not depend on the heavy-particle driving 
forces. Therefore, the first-order electron perturbation function is decoupled from 
the heavy-particles. 

The existence and uniqueness of a solution to eq. (|5.2p is given in the following 
proposition. 



Proposition 5.1. The scalar function 4k given by 



(5.5) ik = -p 5R[mV? c(1) + (M 1 - + iM & )tp, 



Dc(2) 



■d,, 



- M"^ (1) + (M ± + iM )<^ (2) d a 

is the solution to eq. (|5.1j) under the constraints (|4. 1 0[) . where the vectorial functions 
ip^"^, cp^ > °^ 2 \ ifc"^, and ifc"^ are the solutions of the problem 

(5.6) ?c(^ (1) ) - *£, 

(5.7) {J e + i\B\5l°)(^ m )=^ 
where U| e (u) = q e u, under the constraints 

(5.8) «/.V (1) ,^» e = 0. « 6 {1,2}, 
(5-9) UV W .#>>e = 0, ^{1,2}, 



Proof. By linearity and isotropy of the linearized Boltzmann operator 3 e , the de- 
velopment (|5.2p of can be followed through for 4k as well 

1 



4k = -Pc4> c c -d c - </> ° -d x ^— 

The functions cj>%, /i G {-D c , A c }, are now vectorial and satisfy the equations 
(5.10) Jc^c 1 ) + q e C c ^B-dcM = 

and the scalar constraints 

(5-11) <(/ c °0^)) c = 0, I €{1,2}. 

We seek a solution 4>c m the form 

4>» = <^ (1) C e + 4%WC e AB + 4Z {3) C C -B B, 

where 4e^\ <?!>e and 4>e^ are scalar functions of C c -C , (C c -B) 2 and .B£?, since 
4>c must be invariant under a change of coordinates. Substituting this expansion 
in (|5.10p , and using isotropy, eq. (|5.10p splits into three separate coupled equations 

(5.12) 3 C (4Z (1) C C ) - q e B-B4^C c = *£, 

(5.13) %(4>^C e AB) + q e 4^C c AB = 0, 

(5.14) %{4>^C C -B B) + q e C c -B4^B = 0. 

Further simplification is now obtained if, instead of three real quantities 4>e \ 



and 



M3) 



, we introduce one real and one complex unknown defined by 



38 BENJAMIN GRAILLE, THIERRY E. MAGIN*, AND MARC MASSOT 

Upon introducing ip£ (1) = <^ (1) C C and <^ (2) = </^ (2) C c , cqs. (f5~12]) . 1(535]) . and 
(|5.14p can be conveniently rewritten in terms of these new functions 

(EZD (3e + l|B|??)(^ (2) )=^, 

Furthermore, the constraints (|5.11|) are easily rewritten in the form 

m «/eV (1) ,#»e = 0. ^{1,2}, 

(EH) ((/eV^ (2) ,^)) e = 0, ^{1,2}. 

Moreover, expression (|5.5p for </t is immediately obtained using the recombination 
formula 

tf£ = M"^W + M- L $K(^<: 2 )) - M 3(y.^ 2 )). 

□ 



The structure of the integral equation (|5.6[) under the constraints (|5.8[) is clas- 
sical and the structure of equation ()5.7j) under the constraints (15. 9|) is similar in a 
complex framework. More specifically, the operator 3^ + 11-613^° and the associated 
bilinear form a(u, v) = ((u, (CF + i|.B|3^ c )v}} o , defined on the proper Hilbert space 
of complex isotropic squared integrable functions associated with the scalar product 
[•, •], are such that |a(u, u)| > [u, u], which yields existence and uniqueness thanks 
to the constraints. Moreover, from the isotropy of the operator £F e , the expressions 
ipe^ and tpe^ cannot be functions of {C c -B) 2 as shown in (25] . 

We further introduce the electron bracket operators [•, •] and ((•, •)) associated 
with the two operators 3e and For any £ e and £ e , we define 

Ka, Cel e = ((f%, ?e(Ce))) e , ((&, Ce)) e - \B\ (f%, 2?(C0» e - 

These bracket operators develop into 



KcCeL = lE n ^ / <r ej (\C e \ 2 ,u>-e)\C e \f e (\C e \e) 



[e e (|C e | e )-e e (|C e |a;)]©[Ce(|Ce|e)-Ce(|Ce|a;)] du,ded|C e 

Ceel|Ce ~ Cel\fc fcl 



(& + U ~ £e - &)©(& + Ccl - a - &) dwdC e dC e l, 



and 



U*,&\ = \B\q e J / c °e c 0CedC c 



The bracket operator [•,•]„ is hermitian [£ c ,Ce] e = [Cc;£c] c , positive semi-definite 
|^ e ,^e] e > 0, and its kernel is spanned by the collisional invariants, i.e., [£ ,£c] e = 
implies that £ c is a (tensorial) collisional invariant, or in other words, all its 
tensorial components are in the space I c . The bracket operator ((•, ■)) is hermitian 
((Ce,Ce)) e = ((Ce,Ce)) e and negative definite ((£ e ,£e))e < if £ c ¥= 0. 



Remark 5.1. In the limit case B tends to zero, expression (J575J) for the first-order 
electron perturbation function reduces to an isotropic form. We prove indeed that, 
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for /i 6 {£>c,A }, ifc^ does not depend on the magnetic field and that (fe^ 
converges to ipc^ for a vanishing magnetic field. 

The expression of the electron diffusion velocity is given in the following propo- 
sition. 

Proposition 5.2. The electron diffusion velocity V c reads 

(5.15) V c = - (Dl4 + Did^ + D%d°) 

- (6f{d x \nT c )W +et L {d x \nT e ) 1 - + 9 C ® (d x \nT c ) Q ) , 
where the diffusion coefficients and thermal diffusion coefficients are given by 

dL= i Pc t c am^ (1 \<^ (1) L $ =-lM h y^ x \^\, 

(5.16) d± = l Pc T c M h l^\^\, e=-|MU^ (2) ,^ (2) ] e , 



3 t 



Note that the previous expressions are real, in particular for 9^ and 9° Q , although 
functions cp^ e ^ and y>c°' 2 ' are complex. 

Proof. Using definition (|4.16[) of the diffusion velocity V c and expression (|5.4[) of 
*f c yields 

V e =T e M h ((*°%f e %)) e . 
Further substituting expansion (|5 . 5|) into the latter equation, and using isotropy, we 
obtain expression (|5.15j) for the diffusion velocity where the transport coefficients 

are defined by D l L - lp e T e M h ((f e °<p? •))„, = -|M A <(/ V^ (1) , *? e » e , 

Di + iUg = ip c T c Af,((/ c Vf c(2) ,*fi) e , ^+C = -iM,((/ e V^ (2) ,*f c )) c . 
Eqs. ([531) and dUTJ) for n = D c classically yields [25l [27] 



z?i = ip e r e M,[^ (1) ,^ (1) ] , 



3 ^"« Hire i rc j C ' 

e + i0 e es = 4^(i^ (2) , ^ (2) l - 1 ((^- (2) , ^ (2) ) 

As the bracket operators [-, •] and ((■. ■)) are hermitian, we immediately conclude 
for expressions of D e ' e , D^, £>® and 6*c". Concerning 9^ and , we use the 
imaginary part of eq. (|5.7p for /i £ {-D c , A c }, so that 

(5 e + i|B|3*)(^ 2 ))] =0, M e{^ c ,A c }. 

Taking the scalar product of the previous equation with ife , \i £ {-Do, A c } yields 
the four following relations 





-((^ (2) 


,^ (2) 


))e = 


0. 




-((3^ (2) 


,9^ (2) 


))e = 


0, 


fo (2) )Qv5 A s (2) le _ 


-((^ (2) 


,3^ (2) 


))e = 


0, 




-((5ft^ (2) 


,%>^ (2) 


))e = 


0. 
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Then, a direct calculation implies that 

n^ (2 \^ {2 \ = -A e * x > n^ (2 \^ {2 \ - 

^ {2 \v¥ 2 X= 0, 3((^ (2) ,^ (2) )) e = 0, 

so that 9^ = ~\M h {^ {2 \ y^ (2) ] c and 0*> = \M h {{^ {2 \^ {2) )) c . □ 

An alternative form of the diffusion velocity is given by 
(5.17) 

V c = -Dl (4 + ^(fl.lnT.)") - Di (d£ + xf (dJnTe) 1 - + X T(d x lnT e )°) 

- D® (df + X f (O.lnTe)® - xC (9.1nT c )^) , 
where the real thermal diffusion ratios Xe , xS^; Xc are defined by the relations 

(5.18) 9°J = Dlx c c K 9^ + i9°° = (Di + i£®)(xf + i*f )■ 
Then, the viscous tensor is calculated in the following proposition. 
Proposition 5.3. The electron viscous tensor vanishes, i.e., 

(5.19) n c = o. 

Proof. Using definition (|4.37[) of the stress tensor and expression (|5.5p of <fc , one 
readily obtains that II C = 0. □ 

The electron heat flux is given in the following proposition. 

Proposition 5.4. The electron heat flux q e reads 

(5.20) q c = - [At^Te)!! + X^dxTc) 1 - + A®(d x T c )® 

- Pc (e c J4 + e% x di + ef>d°) + Pc h c v c 

where the partial thermal conductivities are given by 

U = ^M^e e(1 We° (1) ] e , 

(5-21) Xi= ^M fc [<^« 

A® =-3^M,((^ (2 \^ (2) )) e . 

Proof. Using definition (|4.30[) of the heat flux q e and expression (|5.4j) of yields 

q e =p e h e V e -M h ((^,f e %)) e . 
Further substituting expansion (|5.5[) into the latter equation, and using isotropy, 
we obtain expression (|5.20[) for the heat flux q c where the transport coefficients 9t" , 
8c 1 - , 9° & are given in eq. (|5.16[) and the partial thermal conductivities Ae , A^ , A® 

are defined by X| = ^M h ((f e °^ {1) , *J*» , X^+iXf = ^M^/eVe^, **•»„. 
Eqs. (pTrJ)) and (pTT)) for [i = A c classically yields [231127] 

Al = 3^M h [^ (1) ,^ (1) ] e , 
Ae + iA® = ^M,([^( 2 ),^( 2 )] e -i((^( 2 ),^( 2 ))) e ). 
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As the bracket operators [•, ■] and ((•, ■)) are hermitian, we immediately conclude 
for the expressions of Xl, X x , and A®. □ 

Using the thermal diffusion ratios defined in eq. (|5.18j) . the electron heat flux is 
rewritten 

(5.22) q e = - [\l(d x T e )U + X^d^T^ + Xf(d x T e )° 

+ Pc (xc 11 ^ 11 + xf % X + X c °\f) + Pch c V c , 
where the thermal conductivities As , A^ , A® are real quantities given by 



All = A 1 ' - — £>" Y e|, Y e| 



T, 



Xc x 



X^ + i\f = + iXf - ^-{Di + i£>® )( X f + ixf )(xf + ixf )• 

The positivity properties associated with the heat flux and the diffusion velocities 
can be written with the help of the mass-energy transport matrices 

(5.23) Al= p*7 I, ] , A ± -'~ A ° ° c » ,0_/- A c <> e 



x -l^cc / \ c ^ee/ \ w c ^ cc 

Proposition 5.5. Considering any two-dimensional real vectors cc", ce\ and x® , 
the two following inequalities are satisfied 

(5.24) (Ajja; 11 ,^ 11 } > 0, 

(5.25) (A^x ± ,x ± ) + (A^x & ,x & ) + (A^ 1 ,^) - (Afx Q , x x ) > 0. 

Proof. Introducing »H = (xf,^)) = (^fSa^j"), and x & = (xf,xf), expressions 
()5 . 16(1 and (|5.2f |1 for transport coefficients yield 

(A ll x ll^ll ) = MkM, [y (i) )2/ (i) L 

(A^x x ,x x ) + {A^x Q ,x & ) + (Afx^^®) - (A * ,^) = B^| y (2) ; y(2)j oi 

with 



ly (i) - Ju,D D (l) i_JI.X(i) 



j/ 2 > = (^ + ^ 2 °)^ (2) - {xi + ix )^ (2) . 

Inequalities (|5.24[) and (|5.25[) are then obtained thanks to the positivity of the 
bracket operator [■, -] c . □ 

Remark 5.2. In the limit case B tends to zero, the behavior of the transport 
coefficients can be investigated. We formally prove that the matrix A\ does not 
depend on the magnetic field, that A^ converges to A\, and that Af vanishes. So 
that we obtain in the limit case the same contribution as with zero magnetic field 
(b = 0) for the electron diffusion velocities and heat flux [55] . 
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5.3. First-order heavy-particle perturbation function. The first-order per- 
turbation function 4h — {<k ) ig jj is solution to eq. (|4.24() . i.e., 

H 

and satisfies the constraints (|4.25|) . where = —$>f(f®)/ 'JJ° ! , i e H. After some 
lengthy calculation based on the expression (I4.8[) of tf, the Euler eqs. (|4.26p . (I4.28[) . 
and (|4.29p . theorem 13. 1[ and coroHarv l3.2[ one obtains 

(5.26) = -*l.d x v h - Ph V *f - *? h -0 x f^) - *f (T a - T h ), 



where 



(5.27) 



- ^(C^Q-iQ-Ql), ieH, 

1 /5f 

2 
77 



= T^(fT h -iTO,;C r Q)a, ZGH, 



Quantity pj = njTft stands for the partial pressure of species i £ H. A linearly 
independent family of diffusion driving forces is also introduced 

(5.28) = — d x pi E F ie , i G H. 

Ph Ph Ph 

The average electron forces acting on the heavy particles belong to the category of 
the diffusion driving forces and allows for a coupling between the heavy particles 
and electrons. Expression of given in eq. (|5.5|) and definition (|4.20p implies that 
F ie , i G H, is proportional to the electron diffusion driving force and the electron 
temperature gradient. Thus, the heavy-particle transport fluxes to be derived are 
also expected to be proportional to the electron forces. 

The existence and uniqueness of a solution to eq. (|5.26p is then established in 
the following proposition. 

Proposition 5.6. The scalar functions family (ffo = ) ieH , given by 

(5.29) <k = -4>l.d x v h - Ph J2 <t>? -d s - tf h d x - 4>f{T e - T h ), i e H, 



is the solution to eq. (|5.26p under the constraints (I4.25P , where the tensorial func- 
tions family = (<f>i) ie n, the vectorial functions families 4> h 1 = (<p i J ) ieH . j G H, 

and (fih — (4 > i h )i £ n? an d t ae scalar functions family </>® = (0f) igH are the solutions 
of the problems 

(5.30) ?,(<#;) = *f , i G H, 

under the scalar constraints 

(5.31) ((//,°C^»fc = 0» iG{l,...,n H + 4} > 
with ,u e {r?, (Dj) je}i ,X h , &}. 
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Proof. By linearity and isotropy of the linearized Boltzmann operator 5i, the de- 
velopment of can be followed through for <\\ as well 

A = -<fi:B m v h - p h Y, <t>? 4- & h d * (jr) ~<t>? (To- T h ), i G H, 



where the functions families 0^, for /j, 6 {ry, (Dj)j e n, \ h , 6} satisfy eq. (|5.30[) under 
the scalar constraints (|5.3ip . We seek a solution in the form 

4>'j = (C4»a - la-ai) , » e h. 

Quantities , /i G {??, (Dj)j E h> an d <^f are scalar functions of Q-Q, for 
i £ H, since 0^, // G {77, (Dj)jgn, A/,, 9} must be invariant under a change of coordi- 
nates. Uniqueness of the solution is readily proved based on the linearity property 
of the operator 3"^, its kernel given in property 14.21 and the constraints (|4.25p 
satisfied by fa. □ 

We further introduce the heavy-particle bracket operator [•,•]/, associated with 
the operator 3^. For any £,h,(h, we define 

[&.Cfc] h = ({//?&, 5fc(Cfc)» h - 

The bracket operator develops into 

[&,a] fc = i E / ^//te+G-d-^)©^ + - a - cpia-QKdwdqdQ. 



The bracket operator [•, -J h is hermitian C^J^ = [On 6il/u positive semi-definite 
l£h, £hjh — 0) an d its kernel is spanned by the collisional invariants, i.e., Oil;, = 
implies that ^ is a (tensorial) collisional invariant, or in other words, that all 
its tensorial components are in the space 2),. The expression of the heavy-particles 
diffusion velocities is given in the following proposition. 

Proposition 5.7. The diffusion velocity of species i G H reads 

(5.32) X = - E D ijdj - e?d x lnT h , 

where the diffusion coefficients and thermal diffusion coefficients are given by 
D l3 = ± Ph T h M h \<p% , 4% ] i,j g H, 



(5.33) 



3 J 

Proof. Using definition (|4.42| of the diffusion velocity and expression (|5 . 2T[) of , 
i G H, yields 

K = rfcM & «*^,/ h %)) h , i G H. 
Further substituting expansion (|5.29p into the latter equation, we obtain expres- 
sion (|5.32p of the diffusion velocities. □ 

In particular, the heavy-particle diffusion velocities are thus proportional to the 
electron driving force and electron temperature gradient through the F ic contribu- 
tion to di, i G H. Kolesnikov [34] has already introduced electron heavy-particle 
diffusion coefficients and thermal diffusion coefficients and ratios to couple the 
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(5.35) 



heavy-particle diffusion velocities to the electron forces. Therefore, we propose 
to refer to this phenomenom as the Kolesnikov effect for the heavy particles. 

^From the properties of the bracket operator, we infer that the diffusion matrix 
D is symmetric. Moreover, an alternative form of the diffusion velocities is given 

by 

(5.34) ^ = -^^(^. + ^111^), zeH, 

where the thermal diffusion ratios are defined from the relations 

Then, we introduce the tensor 

2 

S = [d x v h + [d x v h ) T ] - -d x v h I, 
in order to express the viscous tensor in the following proposition. 
Proposition 5.8. The heavy-particle viscous tensor reads 
(5.36) n h = - Vh S, 

where the shear viscosity is given by 

(5-37) ^ = ^lM h - 

Proof. Using definition (|4.43|) of the viscous tensor and expression (|5.27|) of 
yields 

n h = T h {(*l,f°ct> h )) h . 

Further substituting expansion (|5.29[1 into the latter equation, we obtain expres- 
sion (|5.36[) of the viscous tensor. □ 

The expression of the heavy-particle heat flux is given in the following proposi- 
tion. 

Proposition 5.9. The heavy-particle heat flux reads 

(5.38) q h = -X h d x T h -p h J2 9 j<*3 + Pi h i V v 

where the partial thermal conductivity is given by 
(5-39) h = ^M h \<j>l\4>l h l h - 

Proof. Using definition (|4.45|) of the heavy-particle heat flux and expression (|5.27D 
of *f h yields 

q h = -M h ((*l\fUk)) h + \T h J2 rift. 

Further substituting expansion (|5.29| into the latter equation, we obtain expres- 
sion (|5.38| of the heat flux. □ 
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Using the thermal diffusion ratios defined in eq. ()5.35l) . the heavy-particle heat 
flux is rewritten 

(5.40) q h = -X h d x T h + p h x)V 3 + J2 Pi h i V v 

jGH jGH 

where the thermal conductivity is given by 

(5.41) Xh = X h -n h J2^Xj- 

jen 

The positivity properties associated with the heat flux and the diffusion velocities 
can be written with the help of the mass-energy transport matrix 

Ah = (W Xh 

h [men ('hjh^nj- 

Proposition 5.10. The heavy particles mass-energy transport matrix A h is sym- 
metric, positive semi-definite, and its kernel is onedimensional and spanned by the 
vector [0, (pi)ieB.] T ■ 

Proof. We consider a vector x written in the form x = [xx h , (xi)i£n] T and introduce 
the family y h = (y,) ieH given by 

ieH Ph±h 

Expressions (|5.33|) and (|5.39[) for transport coefficients yield 

(A h x,x) = \p h T h M h {y h ,y h \ h . 

The positivity is then obtained thanks to the positivity of the heavy-particle bracket 
operator [•, ■], . Moreover, using the scalar constraints (|5.3ip . that imply that y h 
is orthogonal to the collisional invariants, the quantity (A h x, x) vanishes if y h is a 
collisional invariant, consequently if y h — 0. Finally, the linear rank of the family 

(4>^ h , (pP 1 , . . . , cjjf'^ 1 ) is exactly rfi because it is the rank of the corresponding right 

member (^f^ h , S&f 1 , . . . , S&f" 11 ). We then conclude that y h — if and only if x lies 
in the space spanned by the vector [0, (pi)i£n] T . □ 

5.4. Second-order electron perturbation function. The second-order pertur- 
bation function </> e 2 is a solution to eq. (|4.33[) . i.e., 

(5.42) J c (4 2 ) + q e dcM)-C^B = * 2 , 
and satisfies the constraints (|4. 34[) . where 

*c = (-4°(/ e : <k ) + M/c°4 , f e % ) + 3 C °) • 

Introducing second-order heavy-particle diffusion driving forces df — —n^, i G H, 
one obtains after some lengthy calculation 

where ^ is a scalar function of C -C c , and 

r^* = ±{c c ®c c - §c e -qj), 



(5.43) 



= ^Q^(\C C \ 2 )\C C \C C , i e H. 
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The coupling of the electrons with the heavy particles occurs in the integral equation 
for the second-order perturbation function through the df forces, i € H. Thus, 
the second-order electron transport fluxes to be derived are also expected to be 
proportional to the heavy-particle forces. 

The complete resolution of eq. (|5.42|) is not necessary since we only need to 
express the second-order transport fluxes V c and q 2 in terms of bracket operators. 
Consequently, we only have to examine the contribution of the two vectorial terms 
and *f«, i £ H. 

Proposition 5.11. The scalar function given by 

(5.44) $ = -e:9A-p e ^»[M«^W + (M x +iM )^( 2 ! 

ieH 



•d? 



is the solution to eq. (|5.42[) under the constraints (|4.34[) . The vectorial functions 
ipe > fe\ i £ H, are the solutions of the problems 



(5.45) WPe iW ) = ^e i , 

(5.46) fa + i\B\3?) =9?, 

under the constraints 

(5.47) <(/ e Vf 4(1) ,^» e = 0, I6{M}, 

(5.48) ((/eV <(2) ,^)) e = 0, ZG{1,2}. 
TTie tensorial function (f>2° verifies 

and the function 4> 2 is a scalar function of C c -C c and (C c -B) 2 . 

Proof. The proof of this proposition is identical to the one of proposition 1 5 . 1 1 since 
eqs. (|5.1[) and (|5.42|) for fa. and <4 2 only differ with their second members. □ 

The expressions of the second-order electron diffusion velocity and heat flux and 
of the average electron force are given in the following proposition. 

Proposition 5.12. The second-order electron diffusion velocity V 2 is given by 
(5-49) V e 2 = -J2( D i d ? 1 + Deid 2 ^ + £>Sd£ ). 



ieH 



The second- order electron heat flux q\ reads 

(5.50) q \ = - Pc £ (C"d- 11 + 0? X df + efd^ ) + Pc h c V 2 . 



The average electron force F ic acting on heavy particles i s H is given by 
(5.51) F lc = ~(Dl4 + D$d£ + D°d?) 

~ T^^t(d x \nT c f + O^id^hxT,) 1 - + 8f(d x \nT c 
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iffusion 


coefficients read 




ef = - 


1 Ti/r [T £>>(!) ^b(1)t 

-\M h l(p c y >,(p c '] c , 


i e H, 


e^ = - 


1 , r [T £>j(2) A e (2) T 


i e H, 


ef = 


1 n r V Di(2) K(2)\\ 


i e h. 


real, in 


particular for D^-, D®, 


#f - 1 - , and 



(5.52) 

4= i Pc T c A/,[^ (1) ,^ (1) ] c , 
DS = -iATeM h ^ W ,^ (a) 3) e , 



!6H, although functions cpe e , <p^ an d V?* a^e complex. 

Proof. Using definition (|4.38|) (respectively (|4.55|) and (|4.20p ) of the second-order 
diffusion velocity V 2 (respectively the second-order electron heat flux q 2 and av- 
erage electron force F ie , i £ H), the same proof as that of proposition 15.21 yields to 
conclude. □ 

Remark 5.3. The term <$!, e :d SB Vh of eq. (|5.44[) contributes to a second-order elec- 
tron momentum relation not investigated here. 

Remark 5.4. To the authors's knowledge, it is the first time that such second- 
order transport coefficients are rigorously derived from a multiscale analysis. The 
second-order electron diffusion velocity and heat flux are thus proportional to the 
heavy-particle diffusion velocities. That is the Kolesnikov effect for the electrons. 
However, it is important to mention that the second-order electron transport fluxes 
shall not be confused with the Burnett transport fluxes appearing in second-order 
macroscopic equations |25) . 

6. Conservation equations 

We review the heavy-particle Navier-Stokes eqs. (|4.52l) - (|4.54p and electron drift- 
diffusion eqs. (|4.60[) and (|4.61[) . We also derive a total energy equation and an 
entropy equation. Then, we introduce a conservative formulation of the system of 
equations. 



6.1. Mass. The species mass conservation equations read 
(6-1) d t p c + dJp c {v h + 1 ^(V c +eV c 2 )) 



0, 



(6.2) 8tPi + d w -[p i (v h + 1 faV i )\=0, ieH. 

Summing eq. (|6.2p over ieH and using the constraint SjeHft^j = given in 
eq. (|4.25|) . a heavy-particle mass conservation equation is obtained 

(6.3) d t ph + d x -(p h v h ) = 0. 

The heavy-particle mass is conserved in the mean heavy-particle velocity referential. 
Then, adding the electron drift eq. (|6.1[) to eq. (|6.3p and using eq. (|2.27p . i.e., 

pv = p h v h + e 2 pe(v h + j^(V c + eV 2 )), 

a conservation equation of global mass p — ph + £ 2 p e is also established 

(6.4) d t p + d x ipv) = Q. 

The global mass is conserved in the hydrodynamic referential, although the trans- 
port fluxes are calculated in the mean heavy-particle velocity referential. It is the 
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only place where the difference between the global hydrodynamic velocity and the 
mean heavy-particle velocity of the order of e 2 , plays an essential role. It is another 
evidence of the coherence of our formalism compared to other approaches found in 
the literature. 

6.2. Momentum. The momentum conservation is expressed by 

(6.5) d t (p h v h ) + d x -{p h v h ®v h + 1 ^pl) = -j^dx-Uh + j^nqE+lSboIo + Sbill^B. 

A flow kinetic energy is obtained by projecting the previous equation onto the mean 
heavy-particle velocity 

(6.6) d t {\ph\vh\ 2 ) + d x - v h {±p h \v h \ 2 + ^p) = -jpp d x v h - ■j^v h -d a> -Tl h 

+ -^pnqE-Vu + v h -(5 bo l + S bl I)/\B. 

h 

6.3. Energy. The electron energy equation reads 

(6.7) d t {p c e c ) + d x - {p c e c v h ) = -p c d x -v h - jj^d x - (q c + eq 2 c ) 

+ ^ (J e + eJ 2 ) -E' + 5 bo eM h J e -v h AB + AE° C + aAE], 
and the heavy-particle energy equation reads 

(6.8) dt(pheh) + d x -(p h e h v h ) = -(p h I + eU h ):d x v h - j^d x -q h + j^J h -E' 

+ AE° h +eAEl 

So that a global energy equation is derived by summing eqs. (|6.7[) and 

(6.9) d t (pe) + d x - (pev h ) = -(pI + eU h ) 

■ eJ 2 + £j h ) -E' + S b0 £M h J e -v h /\B, 



■d x v h - jtf-d x -Q 



where quantity Q = q c + eq 2 + eq h is the total heat flux. Finally a total energy 
equation is derived by adding eq. (|6.6p 

(6.io) d t {£) + d x - (Hv h ) = -sd x -(n h -v h ) - -ir h d x -Q + ie, 

where quantity £ = pe + M 2 ph\\vh\ 2 stands for the total energy and H = p£ + p, 
the total enthalpy. The term I E of eq. (|6.10p represents the power developed by 
the electromagnetic field. It has the form prescribed by Poynting's theorem. Hence, 
the first principle of thermodynamics is satisfied. 

6.4. Electron and heavy-particle entropy equations. In addition to the ther- 
mal energy, we introduce other relevant thermodynamic functions. First, the species 
Gibbs free energy is defined by the relations 



(6.11) p e g e = n e T e ln 



Pi9i = n t T h In 



To' Q° e/ 

where the translational partition functions read 
^m°k B T°\ 3/2 n o_ 



UiTl 



II 



( mi T h f /2 Q° h 



i G H, 



(6.12) Q 

Then, the species enthalpy is given by 
(6.13) 



27rm£k B T 



3/2 
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Finally, the species entropy is introduced as 

h c - g, 



(6.14) 



To 



si = — , % e H. 

J-h 



Therefore, the mixture entropy reads ps = X^gs Pj a j' The thermodynamic func- 
tions exhibit a wider range of validity than in classical thermodynamics, introduced 
for stationary homogeneous equilibrium states [26j . Indeed, they are interpreted 
in the framework of kinetic theory by establishing a relation between the thermo- 
dynamic entropy and the kinetic entropy. The latter quantity is based upon the 
distribution functions 



(6.15) S km = f fi \ 1 - In 



(27r) 3 / 2 n° 



/e l - hi 



(2^) 3 / 2 n° 



fc 



da 



Proposition 6.1. The kinetic entropy and the thermodynamic entropy are asymp- 
totically equal at order e 2 , i.e., 



(6.16) 



S 



kin 



ps + 0(e 2 ), 



provided that the distribution functions follow the Enskog expansion given in eqs. 
(1431) and El. 



Proof. Using definition (|6 . 1 5[) and expansions (|4. 1[) and (14. 2|) . the kinetic entropy 
is found to be 



(2^) 3 / 2 71° 



h 



(27r) 3 / 2 n° „' 



h 



dCj + //e°{l-ln 
dCj + e J f e °<k In 



da 



da 



+ 0{e 2 ). 



The first-order term vanishes seeing the constraints (|4.10|) and (|4.25|) . Then, using 
expressions (|4.7[) and (|4.8p and definition (|6.14p . eq. (|6.16l) is readily obtained. □ 



Consequently, a first-order conservation equation of thermodynamic entropy can 
be used instead of a conservation equation of kinetic entropy to ensure that the 
second principle of thermodynamics is satisfied. First, we introduce the heavy- 
particle entropy phSh — XjeH Pi s j ana * derive the entropy equations. 

Proposition 6.2. The electron and heavy-particle entropy equations associated 
with the macroscopic conservation equations (|6.1[) - (|6.8p read 



(6.17) d t ( Pc s c ) + d x - ( Pe s e v h ) + d x iJ?+eJ^) = T C °+£T C \ 

(6.18) d t {phSh) + d x - {p h s h Vh) + sdvJh = T° + eT^, 
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where the electron and heavy-particle entropy fluxes are given by 

(6.19) £ = -L-fa- M jr e ), J c ' = J ^ F (q 2 c ~ Pc g c V c 2 ), 

M h T e M h T c 

(6.20) j£ = j^jT i^lh ~ X] Pi g i V j) ' 

and the electron and heavy-particle entropy production rates by 



-'I ._ 1 A I?0 ^ c A \r L 

1 



(6.21) T c ° = -A£ e u - jj^rdc-K - —d^T^ - Po h e V e ), 



(6.22) T,! = ^AEl - - -^-d x \nT e -(q 2 - p c h c V e 2 ), 

(6.23) T° = ±AEl 

J-h 



Tl = ±AEl - ±Il h :d x v h - J2 -{9 xPi - n i q i E , yV i 

T h T h M h T h f-J p h 

I 1 -tl 

(6-24) - 1 ±^d x \nT h - (q h - £ p^V, 



Proof. Based on the relations 

/ 9c \ 3n c [ 9i\ 3ri.; 

and definition (|6.14|) . one obtains 

9t(p e s e ) + 9 X - (p c s c v h ) = — [d t {p c e c ) + d x - {p c e c v h )] + n c d x -v h 

J c 

- [<9 t p c + d x - (p e « ft )] =f , 
^ c 

dt(PhSh) + d x - (phShVh) = 7TT [dt(pheh) + d x - (phehVh)] + n h d x -v h 

J-h 

-yz\ d tP] + 9 X - (pjv h )}^-. 

Then, using eqs. (|67Lj) , (jOj) . (|6J|) . (f6T8|) . and the relations 

we readily obtain eqs. (|6.17p and (|6 . 1 8[> . with the entropy fluxes given in eqs. (|6 . 19[) 
and (|6.20[) and the entropy production rates given in eqs. (|6.21I) - (I6.24[) . □ 

Adding eqs. (|6.17p and (|6 . 18[) . a global entropy equation is found 

(6.25) d t (ps) + d x -(psv h ) + d x -J = T, 
where the global entropy flux is given by 

(6.26) J = Jl+ejl + eJl, 
and the global entropy production rate by 

(6.27) T = T c ° + eTl + T ft ° + eT^ 1 . 
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Proposition 6.3. Defining Xh = {d^bxThi d\, ■ ■ ■ , <2„h) t , x\ ~ {{d x \nT c )W , d\) T , 
and = ((c? x lnTc)^, d c ) T , where 

~ d * = ** ~ n ^( D % d ? + ^ e0 (^lnT e )®), 

Ph J- c 

the global entropy production rate T defined in eq. (|6.27[) can be rewritten in the 
following form 



where the matrix A ch 



->O/)c0 n ©n© 



= " c -e— — V D 



is a perturbation of the mass-energy transport matrix A^ defined in eq. (|5.23p . In 
particular, the global entropy production rate is nonnegative provided that e is small 
enough and the collision frequencies Vi e , i S H, are nonnegative. 

Proof. Expression (|6.28p is obtained after some lengthly calculation based on the 
expressions of the diffusion velocities V c , % 2 , \£, i G H, heat fluxes q e , q%, q h , vis- 
cous stress tensor Uh, energy exchange terms AE®, AE*, AE®, AE\, and average 
forces F ie , i 6 H, given in Section [5l 

The positivity of the collision frequencies Ui e , i £ H, (respectively the viscosity 
rjh) immediately yields the positivity of the first term (T c — T^) 2 / (T e Th) Sjgh n j v jc/ m j 
(respectively the second term rjhS:S). Moreover, propositions 15.51 and 15.101 ensure 
that both the following terms jJ h T {A h Xh,Xh) and {A\x\ , x\) are nonneg- 
ative. Finally, introducing y = d^®ip®"^ — (cJ^lnTe)- 1 -®^ ^, the last term is 
expanded as 

(6-29) (A&a£,x±) = [y,y] e -e[z,z] h , 

with 

z t = \ Pc T h M h J2 nj ((y, <p?> (2) ))M? 3 > * £ H. 

We conclude after noticing that the classical term [y, y\ is nonnegative and van- 
ishes if and only if y = thanks to the scalar constraints (|5.8[) and (|5.9p . □ 

In the whole general case, we are not able to write the entropy production rate 
as a sum of nonnegative contributions independently of the value of e. However, it 
is the case for vanishing magnetic field, the last term (A^ h x^, x^r) given in (|6.29|) 
being nonnegative as soon as the magnetic field disappears — the matrix A^ h reduces 
then to A e — and that for any value of e. 

The nonnegativity of the global entropy production rate implies that the second 
principle of thermodynamics is satisfied. This statement could be equivalently 
formulated by means of a H-Theorcm. Besides, the electron and heavy-particle 
temperatures must be equal when an equilibrium state is reached. Provided that 
the collision frequency i/ ic , j e H, is positive, the quasi-equilibrium states described 
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by the Maxwell-Boltzmann distribution functions given in eqs. (|4.7|) and (14. 8| create 
some nonnegative entropy expressed by the term (T c — Th) 2 /(T c Th) X^sh n 3 v ]c/ m ]- 
The latter term vanishes when the electron and heavy-particle temperatures are 
identical. 

6.5. Plasma magnetization. We recall that the intensity of the magnetic field is 
expressed by means of the b parameter used to define the scaling of the nondimen- 
sional electron Hall parameter q° B a t®/m® = e 1 . Three categories of plasmas are 
reviewed in Tabled A value of b < corresponds to unmagnetized plasmas, 6 = 0, 
weakly magnetized plasmas, and 6=1, strongly magnetized plasmas. 





Table 3. 


Magnetic field influence. 


b 


Conservation equations 


Transport properties 


< 









Bulk magnetic force 
Electron magnetic force 


Electron bulk magnetic driving force 


1 


Bulk magnetic force Electron bulk magnetic driving force 
Electron magnetic force Heavy-particle bulk magnetic driving forces 
Heavy-particle magnetic force Anisotropic electron transport properties 



6.6. Mathematical structure. The system of mass, momentum, total energy, 
and entropy eqs. (|6.1[) , (|6.2p . (|6.5[) . (|6.10[) , and (|6.27[) is conservative from a fluid 
standpoint in the variables 

U = [p e , (pi)ieH, PhVh, £■, ps] T ', 

that reads 

(6.30) d t U + d x -F + d x -J-= n, 

with the convective fluxes 

F = [p e Vh, (pi)ienVh, PhVh®v h + — ^pl, Hv h , psv h ] T , 

M h 



the diffusive fluxes 
and the source terms 



n={0,0,^E+(S bo I + 5 bl I)AB, IE, T] T . 

M h 

Then, we extract a purely convective system from eq. (|6.30[) 
(6.31) d t U + d x F = 0', 

where the convective source terms are given by 

M h 
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with the current I' = nqvh, and the entropy production rate 

T c T h m i - 

The purely convective system given in eq. (]6.31[) is rewritten in a quasi-linear form 
(6.32) d t W + A-d x W = n' w , 

by means of the variables 

W=\fi e , (pi)ieH, Wfc, Pe, Ph] T , 

the source terms 

lvl hPh Ph 

and the Jacobian matrices 
(6.33) 








T 



























Vhvl 












IP^u 







\o 





IPhel 








i/G {1,2,3}, 



where the specific heat ratio reads 7 = 5/3 and symbol e v stands for the unit vector 
in the v direction. 

For any direction defined by the unit vector n, the matrix n-A is shown to be 
diagonalizable with real eigenvalues and a complete set of eigenvectors. There are 
two nonlinear acoustic fields with the eigenvalues v^-n ± c, where the sound speed 
is given by c 2 = p/(pf l M?), and linearly degenerate fields with the eigenvalue Vh-n 
of multiplicity + 3. Thus, the macroscopic system of conservation equations 
derived from kinetic theory in the proposed mixed hyperbolic-parabolic scaling has 
a hyperbolic structure, as far as the convective part of the system is concerned. 
Such a property is far from being obvious since the obtained sound speed involves 
the electron pressure and seeing that the rigourous derivation of the momentum 
equation of the heavy particles involves many ingredients throughout the paper. 

7. Conclusions 

In the present contribution, we have derived from kinetic theory a unified fluid 
model for multicomponent plasmas by accounting for the electromagnetic field in- 
fluence, neglecting the particle internal energy and the reactive collisions. Given 
the strong disparity of mass between the electrons and heavy particles, such as 
molecules, atoms, and ions, we have conducted a dimensional analysis of the Boltz- 
mann equation following Petit and Darrozes |46j and introduced a scaling based 
on the £ paramter, or square root of the ratio of the electron mass to a charac- 
teristic heavy-particle mass. The multiscale analysis occurs at three levels: in the 
kinetic equations, the collisional invariants, and the collision operators. The Boltz- 
mann equation has been expressed in the mean heavy-particle velocity referential 
to allow for the resolubility of the first- and second-order electron perturbation 
function equations, as opposed to the inertial referential chosen by Degond and 
Lucquin [19l [20] . Then, the resolubility of the electron and heavy-particle pertur- 
bation functions has been classically based on the identification of the kernel of the 
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linearized collision operators and space of scalar collisional invariants of both types 
of species. The system has been examined at successive orders of approximation 
by means of a generalized Chapman-Enskog method. The micro- and macroscopic 
equations derived at each order are reviewed in Table [2] Depending on the type 
of species, the quasi-equilibrium solutions are Maxwell-Boltzmann velocity distri- 
bution functions at the electron temperature or the heavy-particule temperature, 
therefore, allowing for thermal nonequilibrium to occur. At order e , the set of 
macroscopic conservation equations of mass, momentum, and energy comprises 
multicomponent Navier-Stokes equations for the heavy particles, which follow a 
hyperbolic scaling, and hrst-order drift-diffusion equations for the electrons, which 
follow a parabolic scaling. The expressions of the transport fluxes have also been 
derived: first- and second-order diffusion velocity and heat flux for the electrons, 
and hrst-order diffusion velocities, heat flux, and viscous tensor for the heavy par- 
ticles. The transport coefficients have been written in terms of bracket operators; 
both electron and heavy-particle transport coefficients exhibit anisotropy, provided 
that the magnetic held is strong. We have also proposed a complete description 
of the Kolesnikov effect, i.e., the crossed contributions to the mass and energy 
transport fluxes coupling the electrons and heavy particles. This effect, appearing 
in multicomponent plasmas, is essential to obtain a positive entropy production. 
Besides, it contains, as degenerate case, the single heavy-species plasmas consid- 
ered by Degond and Lucquin for which the Kolesnikov effect is not present. The 
properties of electron and heavy-particle mass-energy transport matrices have been 
established by using the mathematical structure of the bracket operators. In partic- 
ular, the properties of symmetry and positivity implies that the second principle of 
thermodynamics is satisfied, as shown by deriving an entropy equation. Moreover, 
the first principle of thermodynamic was also verified by deriving a total energy 
equation. Finally, the system of equations was found to be conservative and the 
purely convective system hyperbolic, thus leading to a well defined structure. 

The proposed formalism remains valid for collision operators of Fokker- Planck- 
Landau type. These operators can be used to model the charged particle interac- 
tion, instead of Bolztmann operators associated with a Coulomb potential screened 
at the Debye distance. Besides, the explicit expression of the diffusion coeffi- 
cients, thermal diffusion coefficients, viscosity, and partial thermal conductivities 
can be obtained by means of a variational procedure to solve the integral equa- 
tions (Galerkin spectral method [14 ). The expressions of the thermal conductivity, 
thermal diffusion ratios, and Stefan-Maxwell equations for the diffusion velocities 
can be derived by means of a Goldstein expansion of the perturbation function, as 
proposed by Kolesnikov and Tirskiy 35J. Finally, the mathematical structure of 
the transport matrices obtained by the variational procedure can readily be used to 
build efficient transport algorithms, as already shown by Ern and Giovangigli |24] 
for neutral gases, or Magin and Degrez 00] for unmagnetized plasmas. 
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